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Abstract. We present cosmological results from the measurement of baryon acoustic oscilla- 
tions (BAO) in galaxy, quasar and Lyman-a forest tracers from the first year of observations 
from the Dark Energy Spectroscopic Instrument (DESI), to be released in the DESI Data 
Release 1. DESI BAO provide robust measurements of the transverse comoving distance and 
Hubble rate, or their combination, relative to the sound horizon, in seven redshift bins from 
over 6 million extragalactic objects in the redshift range 0.1 < z « 4.2. To mitigate confir- 
mation bias, a blind analysis was implemented to measure the BAO scales. DESI BAO data 
alone are consistent with the standard flat ACDM cosmological model with a matter density 
Qam = 0.295 + 0.015. Paired with a baryon density prior from Big Bang Nucleosynthesis and 
the robustly measured acoustic angular scale from the cosmic microwave background (CMB), 
DESI requires Ho = (68.52 + 0.62) kms”! Mpc^!. In conjunction with CMB anisotropies 
from Planck and CMB lensing data from Planck and ACT, we find Q4, = 0.307 + 0.005 
and Ho = (67.97 + 0.38) kms”! Mpc~!. Extending the baseline model with a constant dark 
energy equation of state parameter w, DESI BAO alone require w = —0.99*0-15. In models 


with a time-varying dark energy equation of state parametrized by w and w,, combinations 
of DESI with CMB or with type Ia supernovae (SN Ia) individually prefer wọ > —1 and 
Wa < 0. This preference is 2.60 for the DESI+CMB combination, and persists or grows 
when SN Ia are added in, giving results discrepant with the ACDM model at the 2.50, 3.50 
or 3.90 levels for the addition of the Pantheon+, Union3, or DES-SN5YR supernova datasets 
respectively. For the flat ACDM model with the sum of neutrino mass ` m, free, combining 
the DESI and CMB data yields an upper limit $m, < 0.072 (0.113) eV at 95% confidence 
for a >) m, > 0 (X>) m, > 0.059) eV prior. These neutrino-mass constraints are substantially 
relaxed if the background dynamics are allowed to deviate from flat ACDM. 
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1 Introduction 


The discovery of the accelerated expansion of the universe [1, 2] demonstrated that the 
dynamics of the universe are presently dominated by dark energy, a component with nega- 
tive pressure [for reviews, see, e.g., 3, 4]. Over the past quarter-century, a wide variety of 
cosmological measurements have lent further support for what has become the standard cos- 
mological model: a spatially flat universe with an energy budget today composed of about 
5% baryonic matter, 25% cold dark matter (CDM), 70% dark energy in the form of Ein- 
stein's cosmological constant (A) and smaller contributions provided by massive neutrinos 
and radiation. 

Increasingly accurate data from type la supernovae (SN la) have confirmed and sig- 
nificantly strengthened the original ground-breaking results for the accelerated expansion 
of the universe (e.g., [5-9]). Meanwhile, mapping out the anisotropies in the cosmic mi- 
crowave background (CMB) radiation, starting with the ground-breaking COBE results [10], 
and continuing with increasingly precise measurements [11-13] that led to the full-sky maps 
by WMAP [14] and Planck [15] experiments, as well as mapping of the CMB anisotropy on 
smaller angular scales [16, 17], revolutionised the field of cosmology by providing percent-level 
constraints on key cosmological parameters. Recent cosmological constraints from probes 
of the large-scale structure [18-21], CMB [22, 23], and distance measurements from SN la 
[24-26] have largely confirmed and sharpened constraints on the cosmological model while es- 
tablishing an increasingly sophisticated methodology framework focused on stringent control 
of systematic errors. 

Current constraints on key cosmological parameters are largely consistent with the 
ACDM cosmological model with cold dark matter and dark energy described by the cos- 
mological constant A. However, dark energy dynamics has not been stringently tested. In 
addition, within the ACDM model, tensions have appeared at various degrees of statistical 
significance (the two most popular being referred to as the “Hubble tension” and the “og 
tension"). Such tensions, if not due to unaccounted systematics, may indicate new physics 
beyond ACDM. 

The use of galaxies as tracers of large-scale structure has traditionally played an essential 
role in cosmology and provides important complementary information to that from SN Ia and 
CMB. Galaxy clustering was pioneered half a century ago and further developed and applied 
to larger and better galaxy catalogs in the intervening years [27-37]. It has established itself 
as a key cosmological measurement that provides a direct link between observations and 
properties of dark matter and dark energy. 

'The clustering of matter encodes a preferred scale, the sound horizon at the baryon 
drag epoch of the early universe [38]. This feature, which is imprinted on the matter distri- 
bution of the early universe by physics around recombination and earlier, is stretched with 
the expansion of the universe, appearing at a comoving galaxy separation of rq ~ 150 Mpc. 
Hence rg is a standard ruler, whose length is dictated by early-time physics. In particular, 
the length of the standard ruler may be calibrated with high accuracy, e.g., by CMB ob- 
servations [15]. Since galaxies trace the matter content of the universe, the BAO feature is 
transferred into galaxy clustering, where it manifests as a single localised peak in the galaxy 
correlation function and an oscillatory signature, or ^wiggles", in the galaxy power spectrum. 
Furthermore, the signature is also visible in other tracers of mass such as fluctuations in the 
Lyman-o forest — spectral features that indicate the radial distribution of neutral hydrogen 
clouds between the observer and distant quasars. Because the BAO feature has a distinct 


signature and resides in the linear-clustering regime, measurements of this scale are relatively 
free from systematic errors associated with nonlinear physics. BAO provides a key cosmolog- 
ical probe sensitive to the cosmic expansion history, with well-controlled systematics. Using 
BAO we can test for dark energy dynamics and spatial curvature, and in combination with 
other probes constrain the Hubble constant, the sum of neutrino masses, and the number of 
light species. Several reviews of the BAO as a cosmological probe are available [e.g., 39-42]. 

The accuracy with which the BAO feature may be measured from a galaxy survey is 
principally limited by sample variance and Poisson noise in the galaxy clustering measure- 
ments, necessitating galaxy surveys with effective volume of at least 1 h73 Gpc? [43-45].! 
The BAO feature was first detected in 2005 by the Sloan Digital Sky Survey (SDSS) [49] and 
the Anglo-Australian Telescope Two-degree Field Galaxy Redshift Survey [50]. Subsequent 
measurements, leading to distance determinations accurate to within a few percent, were 
performed using the SDSS-III Luminous Red Galaxy Sample [51, 52], the WiggleZ Dark En- 
ergy Survey [53-55] and the 6-degree Field Galaxy Survey [56, 57]. Further extensions of the 
SDSS yielded more accurate percent-level BAO measurements using the Baryon Oscillation 
Spectroscopic Survey [BOSS, 58-60] at z < 0.7 and its extension [eBOSS, 61-63] at z > 0.7. 
BAOs have also been detected within the Lyman-o forest mapped by quasar surveys, allow- 
ing distance and expansion measurements at z > 2 [48, 64-67]. Prior to DESI these different 
BAO measurements, considered together, had mapped out the cosmological distance-redshift 
relation with 196-296 accuracy at a series of redshifts in the range z < 2.5 [42, 68]. 

The Dark Energy Spectroscopic Instrument (DESI) is carrying out a Stage IV survey [69, 
70] that was designed to significantly improve cosmological constraints through measurements 
of the clustering of galaxies, quasars, and the Lyman-o forest. DESI is conducting a five-year 
survey over 14,200 square degrees in the redshift range 0.1 < z « 4.2 with a spectroscopic 
sample size that will be ten times that of the previous SDSS surveys. DESI covers six 
different classes of tracers, including low redshift galaxies of the bright galaxy survey (BGS), 
luminous red galaxies (LRG), emission line galaxies (ELG), quasars as direct tracers, and 
Lyman-a (Lya) forest quasars to trace the distribution of neutral hydrogen. Additionally, 
a sample of stellar targets is being observed to a high density in an overlapping Milky Way 
Survey [71] to explore the stellar evolution and kinematics of the Milky Way. For cosmology, 
DESI is designed to impose stringent constraints on both the expansion history and the 
growth rate of large scale structure. Promising detections of the BAO feature at the few 
percent level [72] from the DESI early data release [73] have confirmed that DESI is on target 
to meet the top-level science requirements on BAO measurements. Specifically, DESI will 
tightly constrain matter density in the universe, the equation of state of dark energy, spatial 
curvature, the amplitude of primordial fluctuations, and neutrino mass. It will also sharply 
test modifications to the general theory of relativity proposed to explain the accelerated 
expansion of the universe [74-77]. 

In this paper, we report the constraints derived from the measurements of baryon acous- 
tic oscillations with the first year of data from DESI as part of a larger series of papers. The 
results of this paper build upon the two-point measurements and validation, as detailed in 
[78]; the BAO measurement from galaxies and quasars as summarised in [79]; and the Ly-a 
forest BAO measurements and validation described in [80]. In a follow-up paper, the cluster- 
ing analysis of the first year of DESI data is also performed over a wider range of scales than 
just the BAO feature including the effects of redshift space distortions (RSD) [81]. The Data 


! An exception is the Lya BAO measurement, limited by the number of spectra and their signal-to-noise, 
with only a minor contribution from cosmic variance [46-48]. 


Release 1 (DR1) of the Dark Energy Spectroscopic Instrument [82], comprising spectroscopic 
data obtained in the first year of observations, will be made public at a later stage. 

The paper is structured as follows: Section 2 summarises our data and methodology, 
including a brief overview of the BAO data analysis and the integration of external datasets. 
Section 3 presents the results of the BAO analysis and discusses internal consistencies within 
the DESI data and external consistencies with SDSS data. Section 4 presents our constraints 
on the ACDM model, including a comparison of our results with external data. Section 5 
discusses DESI constraints on extended dark energy models both independently and in com- 
bination with CMB and SN la data. In Section 6, we explore the constraints on the Hubble 
parameter for a variety of cosmological models and dataset combinations. Section 7 sum- 
marises our constraints on the neutrino sector. Lastly, we summarise our findings in Section 8. 


2 Data and methodology 


In this section, we introduce the cosmological quantities relevant for this analysis, review the 
basics of the BAO as a cosmological probe, then describe the DESI data we analyze, the 
external datasets we combine with DESI data, and the analysis methods we adopt. 


2.1 Cosmological background 


In a homogeneous and isotropic cosmology (i.e., with the Friedmann-Lemaítre-Robertson- 
Walker (FLRW) metric), the transverse comoving distance is 


Du(z) = ZI sinh [vine f. HG zu (2.1) 


Here z is redshift, Ho is the Hubble constant with H(z) the Hubble parameter, Qx is the 
curvature density parameter, and c is the speed of light. In ACDM, Ok = 1— Qu — O4 — OR, 
where Om, Og, and QA are the energy densities relative to critical in matter, radiation, and 
the cosmological constant, respectively, all evaluated at the present time. Since sinh is a 
complete function, the expression above holds for positive, negative, or zero Qg. The Hubble 
parameter in ACDM is 


H(z) = HovVQa(1+ 23 + On (1+ z)4 + 9k (1+ z)? +94. (2.2) 


The Hubble parameter also contains the contribution by massive neutrinos given by their 
energy density relative to critical (2,. They are non-relativistic at low redshifts that DESI 
tracers probe and hence contribute to H(z) like matter, but their contribution becomes that 
of a relativistic species at high redshift. 

BAO measurements depend on the sound horizon at the drag epoch rg. This is the 
distance that sound can travel between the Big Bang and the drag epoch which indicates 
the time when the baryons decoupled. The drag epoch occurs slightly later (at za = 1060 
in the standard model) than photon decoupling (z, ~ 1090) simply because there are so few 
baryons relative to the number of photons. The expression of ra is given by 


oo 
Gs(2) 
„ HG) e) 
where c, is the speed of sound which, prior to recombination, is given by 
C 
Glz) = (2.4) 


pa) 


where pg and p, are the baryon and radiation energy densities, respectively. The speed of 
sound evaluates to approximately cs ~ c/v3 well before recombination, but then decreases, 
and finally drops sharply at recombination. For standard early-time physics assumptions, 
the drag-epoch sound horizon can be approximated with [83] 


147.05 ( Um Rd Ng ( Wh pm (2.5) 
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where wm = Qmh? and wp = Q^? are the matter and baryon physical energy densities today, 
respectively, and Neg is the effective number of extra relativistic degrees of freedom. 

'The sound horizon at the drag epoch leaves an imprint in the distribution of matter 
that serves as a cosmological standardised ruler [46, 84-86]. To illustrate the concept, con- 
sider an ensemble of galaxy pairs at a given redshift z: if the separation vectors of the pairs 
are oriented perpendicular to the observer's line-of-sight, a preferred angular separation of 
galaxies A0 may be observed, measuring the comoving distance Dy(z) = ra/ÃA6 to this red- 
shift. With the separation vector oriented parallel to the line-of-sight, a preferred redshift 
separation Az may be observed, measuring a comoving distance interval that for small inter- 
vals gives the Hubble parameter at that redshift, represented in this paper by the equivalent 
distance variable Dy(z) = c/ H(z) = ra/ Az. BAO measurements hence constrain the quan- 
tities Du(z)/ra and Dy(z)/rg. This interpretation holds true for standard assumptions and 
models not too dissimilar from ACDM, given the statistical power of the data [87]. For those 
BAO measurements in certain redshift bins with low signal-to-noise ratio, we instead quote 
constraints on the angle-averaged quantity, Dy(z)/ra, where Dy(z) is the angle-average dis- 
tance that quantifies the average of the distances measured along, and perpendicular to, the 
line of sight to the observer [49]: 


1/3 
Dy(z) = (zDu(2)? Du(2)) ^. 

BAO measurements are performed at a series of redshifts, allowing constraints to be 
obtained on the cosmological parameters governing the distance-redshift relation including 
curvature and dark energy, and the Hubble constant if external information on the absolute 
BAO scale is provided. 


(2.6) 


2.2 DESI BAO data from Data Release 1 


DESI spectroscopic targets are selected from photometric catalogs of the 9th public data 
release of the DESI Legacy Imaging Surveys [88],? drawn from three optical surveys in the 
grz optical bands: DECaLS using the DECam camera [89] (which includes data from the 
Dark Energy Survey (DES) [90]) south of declination 32.375?, and north of this limit the 
Beijing-Arizona Sky Survey (BASS) [91], and the Mosaic z-band Legacy Survey (MzLS) [92]. 
Four different classes of extragalactic targets are defined: the bright galaxy sample (BGS, 
[93]), luminous red galaxies (LRG, [94]), emission line galaxies (ELG [95]), and quasars (QSO 
[96]). 

Spectroscopic observations of these targets are carried out with the DESI instrument [70] 
mounted on the Nicholas U. Mayall Telescope at Kitt Peak National Observatory in Arizona. 
Each observation field is covered by a "tile", consisting in a set of targets located within 
that sky area [97] and assigned to each of the 5000 fibers in the focal plane of the telescope. 


?https://www.legacysurvey.org/dr9/ 


Each fiber is placed at the celestial coordinates of its assigned target by a robotic positioner 
[98, 99] and guide its light from the focal plane to one of the ten DESI spectrographs. 
DESI observing time is dynamically separated into a “bright” time program (when BGS are 
observed) and “dark” time observations (when LRG, ELG, and QSO are targeted) depending 
on observing conditions. Redshift distributions, exposure times, calibration and observation 
procedures were optimised during a period of survey validation [100] that included a visual 
inspection campaign [101, 102]. The DR1 spectroscopic dataset is built from main survey 
operations starting from May 14, 2021 through June 14, 2022, using an observing strategy 
meant to prioritise depth [103], resulting in 2,744 “dark” tiles and 2,275 “bright” tiles. The 
covered tile surface area is of order 7,500 deg?, just over half of the expected final coverage 
of 14,200 deg?. However, the completeness within this area will significantly increase, as 
one can infer given an expected final number of 9,929 dark and 5,676 bright observed tiles 
[103]. The combined effective volume is expected to increase by more than a factor 3 [79]. 
The observed data are processed by the DESI spectroscopic pipeline [104] on a daily basis 
for immediate quality checks. The redshift catalogs used for this analysis and released with 
DESI DR1 are obtained from a spectroscopic reduction run with a fixed pipeline version 
internally denoted as “iron”. 

Large-scale structure catalogs of galaxy and quasar positions suitable for the clustering 
analysis are built from the redshift and parent target catalogs and their two-point function 
measurements are obtained with all DR1 specific details presented in [78]; a technical overview 
of the general process is presented in [105]. The selection function is defined and correction 
weights are designed to compensate for systematic variations due to the effects of imaging 
anisotropies on the input target samples [106, 107], redshift measurement efficiency [108, 109], 
and incompleteness in target assignment [110-112]. Simulations of the DR1 dataset are 
presented in [113]. 

Studies of the Lya forest are based on a quasar catalog with alternative redshift esti- 
mates to minimize possible biases caused by the same Lya absorption [114]. The method to 
extract Lya fluctuations from the quasar spectra is described in [115], including the masking 
of pixels contaminated by Broad Absorption Line features [116] and Damped Lyman-a sys- 
tems [117]. Below, we provide details on each of the target samples studied and the properties 
of the particular samples that are used for this analysis. 


The Bright Galaxy Sample (BGS, 0.1 < z < 0.4): The nominal target selection (BGS 
Bright) relies on r-magnitude cuts tuned to achieve uniform density across the photometric 
samples. Gaia catalogs [118] are additionally used to remove point-like sources. The resulting 
sample has 854 targets per square degree, which are assigned fibers with high priority during 
bright time observations. 

The nominal BGS sample is flux-limited and high density, with strong evolution in 
redshift. Over 5.5 million reliable BGS redshifts were measured in DR1. However, in order 
to produce a more homogeneous sample, a cut was engineered to produce a sample of roughly 
constant comoving number density of 5 x 107* h3 Mpc”? using k-corrected r-band absolute 
magnitudes from [119, 120] and a redshift-dependent correction for evolution (matching the 
sample used for [73]). The BGS number density is similar to that of LRGs at z = 0.4 and 
high enough to make shot-noise a minor contribution to the BAO statistical uncertainty (see 
[78] for more details). The final BGS clustering sample used for the BAO measurement 
comprises 300,017 redshifts in 0.1 < z < 0.4, with an assignment completeness of 61.6%, 
which is expected to increase to ~ 80% in the finalized survey. 


The Luminous Red Galaxy Sample (LRG, 0.4 < z < 0.6 and 0.6 < z < 0.8): The 
LRG target selection [94] uses photometry from the Wide-field Infrared Survey Explorer 
(WISE) [121] to select red objects with a nearly constant comoving number density of 5 x 
10-4 A? Mpc^? over 0.4 < z < 0.8. An additional cut on the z-magnitude measured within 
the aperture of a DESI fiber further selects targets with high S/N spectra. The obtained 
624 deg"? target sample is assigned fibers with priority higher than ELG targets but lower 
than QSO targets. 

'The DESI DR1 LRG clustering sample used for BAO measurements consists of 2, 138, 600 
redshifts in the redshift interval 0.4 < z < 1.1. The DR1 assignment completeness is 
69.296, which is expected to increase to 9096 in the completed survey. The lower redshift 
bound was chosen to separate the sample from BGS, as most low-redshift LRG targets are 
also BGS targets, while the upper bound was designed to ensure a minimum density of 
10-4 A? Mpc^?. This sample is further split for the clustering analysis into 3 disjoint red- 
shift ranges, 0.4 < z < 0.6, 0.6 < z < 0.8, and 0.8 < z < 1.1 (for brevity, we sometimes refer 
to these as LRG1, LRG2 and LRG3, respectively), however we do not use any BAO measure- 
ment from the highest-redshift LRG bin alone. It is instead combined with the overlapping 
lowest-redshift ELG bin, as described in the combined LRG and ELG sample section below. 


The Emission Line Galaxy Sample (ELG, 1.1 < z < 1.6): ELG targets are selected 
with colour cuts in the (g — r) vs (r — z) plane to prioritise objects with [OII] emission in the 
desired redshift range [95] for secure redshift measurements. Low-redshift objects are filtered 
out with a g-magnitude cut and another cut on the g-magnitude within the aperture of a 
DESI fiber further favors targets with high S/N spectra. While a random 10% fraction of 
the ELG sample receives the same fiber assignment priority as LRGs to facilitate small-scale 
clustering measurements in the redshift range where the samples overlap, the remaining 9096 
of ELGs are assigned fibers at a lower priority than LRGs and QSOs. 

The DR1 ELG sample defined for clustering analysis in [78] comprises 2, 432, 022 reliable 
redshifts in the interval 0.8 « z « 1.6. T'he lower redshift bound aligns with that of the high 
redshift LRG sample, while the upper bound rejects objects whose [OIT] doublet falls outside 
of the spectrograph wavelength coverage. In DR1, the fiber completeness is only 35%, which 
should increase to over 6096 in the final dataset. The sample is split into two disjoint redshift 
ranges for clustering analysis, 0.8 < z « 1.1 and 1.1 < z « 1.6 (similarly referred to as 
ELG1 and ELG2 for brevity). We do not use any standalone BAO measurement from the 
low-redshift ELG bin; it is instead combined with the overlapping high-redshift LRG bin. 


The combined LRG and ELG Sample (LRG+ELG, 0.8 < z < 1.1): The high- 
redshift LRG and the low-redshift ELG samples overlap in the range 0.8 < z < 1.1. These 
two samples are concatenated with inverse-variance weights to obtain a combined LRG+ELG 
catalog, following the methodology presented in [122]. The increased combined density is 
expected to improve the reconstruction efficiency, while the inverse-variance weights are de- 
signed to maximize the measurement of the BAO precision. As described in [79], the ob- 
tained combined LRG+ELG BAO measurement is ~ 10% more precise, while being fully 
consistent with the LRG-alone measurement in this redshift range (and with that of ELG 
alone, although this has a larger uncertainty). We therefore use the combined LRG+ELG 
0.8 < z < 1.1 BAO measurement for the cosmological inference. 


The Quasar Sample (QSO, 0.8 < z < 2.1): The QSO target selection relies on identi- 
fying a flux excess in the near-infrared, which is assessed by comparing near-infrared WISE 


W1, W2 magnitudes to Legacy Imaging Surveys optical grz magnitudes [96]. Sources with 
stellar morphology (PSF type) are selected, and a r-magnitude cut is applied to eliminate 
bright stars and faint targets. Quasars are then selected from the r magnitude and 10 colors 
using grzW1W2 with a Random Forest (RF) algorithm, designed to produce a target sam- 
ple of density 310 deg”?. The resulting QSO target sample is assigned fibers with maximum 
priority. 

The DR1 QSO sample used for BAO measurements in 0.8 < z « 2.1 consists of 856, 652 
redshifts. Quasars at z > 2.1 are not used as part of this sample. Because of their maximum 
priority, the assignment completeness of QSO in DR1 is already high, 87.6%. 


The Lyman-a Forest Sample (Lyo, 1.77 < z < 4.16): The highest-redshift BAO mea- 
surement from DESI DR1 is described in [80], and is obtained from a combined analysis of 
correlations of three different datasets. The first dataset consists of the positions (celestial 
coordinates and redshifts) of 709,565 quasars in the redshift range 1.77 < z « 4.16. In 
addition, the Lya forest in the spectra of high-redshift quasars at z > 2.1 constitute an al- 
ternative tracer of matter density fluctuations. We use two Lyo forest datasets, consisting of 
fluctuations in two different rest-frame wavelength regions of the background quasar spectra: 
428, 403 in region A (1040 < A, « 1205 À) and 137,432 in region B (920 < A, < 1020 À).? 
In order to increase the signal-to-noise of the spectra, quasars identified as having z > 2.1 
after a single epoch of observation are prioritised for up to four observations. 

Table 1 summarises the properties of the different samples used in this analysis, as well 
as the BAO measurements obtained from them using the methods described below. Overall, 
this analysis of the DESI Data Release 1 relies on over 6 million unique redshifts, more than 
twice the number of redshifts considered in the final SDSS cosmological analysis [42]. The 
aggregate precision on the BAO isotropic scale is 0.4996 to be compared to 0.6096 for the 
final SDSS measurements. 


2.3 BAO analysis 


We now briefly describe how DESI BAO analyses of DR1 data implement the standard ruler 
measurement. 

BAO analyses usually need to assume a fiducial cosmology although the resulting dis- 
tance and expansion measurements remain independent of these assumptions in a leading- 
order approximation, unless the assumed fiducial model differs significantly from the model 
to be tested, or the truth [123-125]. First, a fiducial cosmological model is adopted to convert 
tracer angular positions and redshifts to comoving coordinates, and the two-point clustering 
pattern is measured as a function of comoving separation assuming this fiducial model. The 
clustering pattern perpendicular and parallel to the line-of-sight changes in a predictable 
way between fiducial and trial cosmologies according to the Alcock-Paczynski (AP) effect 
[126, 127], allowing a wide range of cosmological models to be tested [128]. 

Second, the theoretical calibration of the BAO scale is encoded by a full early-universe 
matter power spectrum generated assuming a fiducial cosmological model [38, 129]. This 
matter power spectrum is used to produce a template tracer clustering model, indexed by 
the sound horizon at baryon drag, which is scaled by two free scaling parameters—along and 
perpendicular to the line-of-sight—during data fitting. In order to measure the location of 
the BAO feature alone, the model also includes “nuisance” parameters for the broad-band 
shape of the tracer power spectrum, which are marginalized over in the analysis [58]. Residual 


“Region B is treated separately from region A as it is also affected by higher order Lyman lines. 


theoretical systematics arising from nonlinear matter growth, tracer bias, and redshift-space 
distortions are reduced by this procedure, and are generally considered to be present at a 
level of ~ 0.1% [130-133]. The modeling and the procedure for performing BAO galaxy 
and quasar fits have been updated for DESI from the broadband polynomial marginalization 
procedure used in previous surveys, according to the theoretical framework laid out in [134]. 

Finally, the amplitude of the BAO feature imprinted in the early universe is diluted 
by the bulk-flow motion of matter from its true locations in the density field, driven by the 
growth of structure, which changes tracer separations by a few h-! Mpc. This effect may 
be partially reversed using density-field reconstruction to estimate these displacements based 
on the gravitational field inferred from the observed tracer distribution [135, 136]. Applying 
these reversed displacements to tracers sharpens the BAO signature and mitigates the small 
(« 0.396) nonlinear shift of the BAO-scale in late-time clustering, allowing for more accurate 
cosmological inferences [137, 138]. The reconstruction algorithm used has been updated 
relative to that used for eBOSS [139] and the specific choice of algorithm and settings were 
validated as described in [140]. 

For the galaxy and quasar BAO data points considered in this paper we make use of 
the post-reconstruction clustering measurements in configuration space presented in [79]. 
Anisotropic BAO measurements are obtained for LRG and ELG samples, but we restrict 
to isotropic fits for BGS and QSO samples due to the lower signal-to-noise ratio of the 
two-point function measurements. The uncertainties in these clustering measurements are 
described by covariance matrices whose construction and validation are outlined in [141- 
143]. Various contributions to the total systematic error budget arising from theoretical 
modeling uncertainties (which are at most 0.1% and 0.2% of the isotropic and anisotropic 
BAO parameters, respectively [134]), uncertainties due to the galaxy-halo connection (of 
< 0.2% [144, 145]) and observational systematic effects (which are negligible [79, 106, 108]) 
have been estimated in a series of supporting papers. These are added together in quadrature 
as summarised in [79] to compute the total error budget for the BAO measurements in the 
BGS, LRG, ELG, and QSO samples. Together, these steps represent the state-of-the-art in 
the measurement of BAO from galaxy tracers and are described in detail in the companion 
DESI paper [79]. 

In addition to measuring the clustering of galaxies and quasars, we also make use of the 
Lya forest. The Lya BAO measurement is presented in [80], and it combines a measurement 
of the auto-correlation of fluctuations in the Lya forest and of its cross-correlation with 
the density of quasars at z > 1.77. The method to measure and model the correlations is 
described in [146], and it is based on previous Lya BAO measurements from eBOSS [67]. In a 
supporting study [147] we validate the analysis using synthetic datasets, simulated using the 
methodology described in [148]. An important improvement in the methodology is that we 
now take into account the small correlation between the measurements of the auto- and the 
cross-correlation. Finally, the (minor) impact of correlated noise in the Lyo auto-correlation 
is studied and characterised in [149]. Besides the validation tests using synthetic datasets, 
in [80] we also present several data splits and a long list of consistency tests that show that 
our Lya BAO measurement are robust to changes in the measurement methodology or in 
the modeling. For instance, we show that adding an ad-hoc smooth component to the model 
with up to 48 extra free parameters does not vary the BAO parameters by more than 0.196. 

The BAO measurements obtained from the various samples used are shown in Table 1. 
As already mentioned, for the BGS and QSO tracers, we only measure the angle-averaged 
Dy/ra quantity, due to the lower signal-to-noise achieved. For all other tracers, we quote 


tracer redshift Nexaéer Zeg Du/ra Dg/ra r or Dy/ra Lu 
(Gpc”) 
BGS 0.1 — 0.4 300,017 | 0.295 — — 7.93+0.15 1.7 
LRG1 0.4 — 0.6 506,905 | 0.510 | 13.62 + 0.25 | 20.98 + 0.61 —0.445 2.6 
LRG2 0.6 — 0.8 771,875 | 0.706 | 16.85 + 0.32 | 20.08 + 0.60 —0.420 4.0 
LRG3+ELG1 0.8 -- 1.1 1,876,164 | 0.930 | 21.712: 0.28 | 17.88 + 0.35 —0.389 6.5 
ELG2 1.1— 1.6 1,415,687 | 1.317 | 27.79 + 0.69 | 13.82 + 0.42 —0.444 2.7 
QSO 0.8—2.1 856,652 | 1.491 — Bê 26.07 + 0.67 1.5 
Lya QSO 1.77 — 4.16 | 709,565 | 2.330 | 39.71+0.94 | 8.52+0.17 —0.477 — 


Table 1. Statistics for the DESI samples used for the DESI DR1 BAO measurements used in this 
paper. For each tracer and redshift range we quote the number of objects (Niracer), the effective 
redshift (zem) and effective volume (Ve). Note that for each sample we measure either both Dy/ra 
and Dg/ra, which are correlated with a coefficient r, or Dy/ra. Redshift bins are non-overlapping, 
except for the shot-noise-dominated measurements that use QSO (both as tracers and for Lya forest). 


results in Dy /rq and Dy/ra, and provide the value of the correlation between them, r. 


2.3. Blinding framework 


An imperative of science with modern cosmology experiments, and especially for a survey 
with the statistical precision that DESI is able to achieve, is to mitigate against the possibility 
of the results being affected by observer confirmation bias. This can manifest itself through 
the attribution of unexpected results, or results which do not match the prior prejudices of 
the observer, to systematic effects which must be corrected for in the data analysis. Thus 
the observer may—consciously or subconsciously—continue to search for and correct real or 
imagined systematic effects in the data only until they agree with some pre-conceived desired 
result, and no further. 

To mitigate against the possibility of such biases unknowingly entering into our analyses, 
we apply a system of blinding our data, to conceal the true cosmological results during the 
period where systematic errors were being investigated and analysis pipelines finalised. This 
blinding was applied in two different ways: 


e For the discrete galaxy and quasar samples, blinding was applied at the catalog level 
[150], according to the procedure described and validated in [151]. The redshifts of 
the tracers in the catalogs were shifted in such a way to change the position of the 
BAO feature as well as the redshift-space distortion (RSD) signal, and the weights 
applied to the tracers were adjusted in order to blind the measurement of primordial 
non-Gaussianities from the large-scale clustering. True redshifts and weights were not 
made available for the analysis before unblinding. 


For the Lya sample, catalog-level blinding was challenging due to sky lines and Galactic 
absorption features located at known observed wavelengths in the spectra. Instead, 
the apparent position of the BAO peak was shifted by directly applying an additive 
component to the measured correlation functions, as described in [80]. 


In both cases, all analyses used the blinded data during initial testing and systematic errors 
were identified and corrected. Only after the analysis pipeline had been finalised, the system- 
atic error budget determined, and a series of strictly prescribed validation tests passed (see 
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[79, 80] for a detailed description) was the data blinding removed and the unblind results 
revealed. For discrete tracers, it was decided prior to unblinding that the final choice of 
covariance matrices would be made after unblinding, and that the LRG and ELG samples 
would be concatenated in 0.8 < z < 1.1 for a combined BAO measurement in this redshift 
range. 

We mention two unplanned updates in the BAO analysis pipeline for discrete trac- 
ers made after unblinding and decided independently of the actual BAO measurements: 
a completeness correction was fixed in the clustering catalogs, and the implementation of 
the BAO theoretical model was corrected to fully match the theoretical framework laid out 
in [134]. Corrections in the theoretical model changed BAO measurements by less than 
0.30 (differences in most redshift bins under 0.10), and updates to the catalogs by at most 
0.60 (variations in most redshift bins being under 0.20). We also report two minor post- 
unblinding changes in the Lya BAO analysis: a bug fix that led to a < 0.1e change in the 
BAO measurement, and an update in a bias parameter that had no effect on the results. 

Two additional points are worth noting. For the discrete galaxy and quasar samples, the 
analyses related to the full shape of the broadband clustering signal (described in [81]) and 
those pertaining to the BAO results presented in [79] and included in this paper had separate 
pipelines and robustness criteria to be satisfied before unblinding, and were not unblinded 
at the same time. Additionally, it was decided that cosmological model inference from the 
BAO and full shape measurements would remain outside of the scope of the blinding, and 
would be performed once BAO and full shape unblinded results are obtained.* 


2.3.2 Summary of DESI DR1 BAO likelihoods 


The results used in this paper are derived from the BAO measurements and fits described 
above applied to each of the five tracer samples across 7 distinct redshift bins covering a total 
redshift range from 0.1 to 4.2. 

The BGS and LRG samples chosen cover disjoint redshift ranges. As a result of the lack 
of overlap, the correlation between the BAO measurements in these samples is small, and we 
assume it to be negligible and do not include any covariance between them in our cosmological 
analysis. In the 0.8 < z < 1.1 redshift bin, where the LRG and ELG tracers overlap, they 
are analysed together using a single multi-tracer approach, described in [79, 122]. The QSO 
sample overlaps with this redshift bin, and with the ELGs in the range 1.1 < z < 1.6. Ina 
general case where two tracer samples overlap in volume, the correlation coefficient between 
two estimates of the power spectrum obtained from these tracers may be estimated as: 


XıXo 
5b f dV 
ae PRE. J X12 (2.7) 
D A 1/2 1/2? : 
V Var Pi) Var( £5) (f dV X1)! ($ dV X2)¥/ 
B PLP. 2 
where X; = (Ee) for i = 1,2, X12 = (= aati] , my and ng are the mean 


tracer number densities of the two samples and n12 is the common sample density, and Py 
and P» are the corresponding typical amplitudes of the power spectra of the two tracers at 
the BAO scale k ~ 0.1 A Mpc^!. For different tracer types, such as the overlapping ELG 
and QSO samples of relevance in the 1.1 « z « 1.6 redshift bin, there are no common tracers 


“Internal checks of the inference pipeline were performed with synthetic data prior to unblinding. 
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so ni2 = 0. This correlation between power spectrum estimates can be roughly regarded 
as an approximation to the correlation between the BAO scale measurements from the two 
overlapping tracer samples, although we regard this only as a guide and caution against using 
this estimate for any precise combination of results from different tracers. 

For the ELG and QSO tracers in 1.1 < z « 1.6, using representative values in Eq. (2.7) 
gives a very small correlation, C « 0.1. T'he correlation between QSOs and LRGs and ELGs 
in the lower redshift 0.8 < z « 1.1 bin is even smaller. We therefore ignore this correlation 
and treat the QSO BAO measurement as effectively independent of those from LRG and 
ELG tracers. We also ignore the correlation between the BAO measurement from the QSO 
sample (at z,g = 1.49) and the Lya BAO measurement (at zeg = 2.33). Even though 
quasars in the range 1.77 « z « 2.1 contribute to both measurements, in the Lya BAO 
measurement they are only used in cross-correlation with Lya absorption at z > 2. Therefore, 
these measurements would only be correlated by cosmic variance, but both measurements 
are dominated by shot-noise. In this redshift range, the number density of the quasars is 
less than 3 x 107? ^? Mpc^? and given the P(k) at the scales most relevant for BAO is 
< 7500 h-? Mpc?, nP < 1 and any correlation between the measurements will be small. 

As described in Section 2.1, the cosmological quantities measured by the BAO analysis 
in each redshift bin are either the two (correlated) distance ratios (Dy /ra, Dy/ra) or just 
a single one, Dy/ra, depending on the signal-to-noise ratio. The recovered posteriors in 
these quantities are in all cases found to be well approximated by Gaussian distributions. 
Systematic uncertainties are added in quadrature to the Gaussian covariance matrices, which 
we use together with the mean vectors for the cosmological inference presented hereafter. 


2.4 External Datasets 


In this paper, we compare cosmological results from our data to those from a number of other 
recent experiments and provide joint constraints from different combinations of datasets. We 
briefly describe these data and likelihoods below. 


2.4.1 Big Bang Nucleosynthesis and Qh? 


In a flat ACDM background model, the BAO distance ladder results determine Qm and Hora. 
From Eq. (2.5), the sound horizon depends on the physical matter and baryon densities, 
Wm = Oy h? and wp = 9p h2. Therefore, assuming standard neutrino content of the universe 
(Neg = 3.044), we require prior knowledge of Qh? to break the degeneracy between Ho and 
Td: 

Big Bang Nucleosynthesis (BBN) theory predicts the abundances of light elements in- 
cluding D and *He in the early universe, and these abundances depend on the baryon to 
photon ratio 719. Thus the observational determination of the primordial deuterium abun- 
dance D/H [152] and the helium fraction Yp [153, 154] can be used to deduce O» h2. 

The resulting constraints on Qh? depend on the details of the theoretical predictions 
and their treatment of underlying nuclear interaction cross-sections, in particular for the 
deuterium burning reactions. For a review of the details we refer interested readers to [155]. 
A recent analysis [156] makes use of the new PRyMordial code [157] to recompute the predic- 
tions while correctly marginalizing over uncertainties in the reaction rates, and reports the 
conservative constraints 

QyA? = 0.02218 + 0.00055 (2.8) 


Vi2 P4 P: 


In this case and in the limit of a thin redshift slice, Eq. (2.7) simplifies to C = MVP am) (Pati fia)” 


as used in [139]. 
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in the standard ACDM model, and 


QyA? = 0.02196 + 0.00063 (2.9) 


when allowing for additional light relics (ACDM-+ Nog model), where in the latter case there 
is also a covariance between Qh? and Neg which we account for. We adopt these values 
as BBN priors in the following sections. Note that these BBN priors are more conservative 
than some previously proposed in the literature (e.g., [158]); this is a direct result of the 
marginalization over reaction rate uncertainties performed in [156]. 


2.4.2 Cosmic Microwave Background 


The power spectra of anisotropies in the cosmic microwave background (CMB) have been 
exquisitely measured by Planck [15, 159]. A fundamental quantity that is determined from 
the oscillations in these power spectra is the measured sky angular scale of acoustic fluctu- 
ations, i.e. 0, = r./Dw(z.), where rą is the comoving sound horizon at recombination and 
Du(z+) is the transverse comoving distance to the redshift of recombination. This angular 
scale is a very sharp feature and can be measured precisely independently of possible sys- 
tematic effects affecting the high-/ part of the CMB power spectrum (in particular the Ar, 
parameter). 
Within the ACDM model Planck reports [15] 


10062" — 1.04109 + 0.00030, (2.10) 


a 0.0396 precision measurement through the combination of temperature and polarization 
data alone. The central value of 02c! is almost entirely independent of the specific cosmo- 
logical model assumed, although the uncertainties can increase significantly in some model 
extensions affecting early-time physics. This is a consequence of the simple geometrical inter- 
pretation: 6, is the BAO scale imprinted in the CMB anisotropies at recombination. DESI, 
on the other hand, measures the same BAO features imprinted in the galaxy distributions 
at lower redshifts. It is therefore natural to combine DESI BAO results with the constraint 
on 0, from Planck, which can be regarded as a minimal and highly robust purely geometric 
version of the full CMB information. To achieve this, we implement a Gaussian external prior 
on the quantity 1000, with mean 1.04110 and variance 0.000537. The mean here matches the 
reported value in [15] when including CMB lensing, and is very close to that in Eq. (2.10). 
We have conservatively increased the width of the prior by ~ 75% relative to the baseline 
result in order to accommodate the increased uncertainties in models beyond standard flat 
ACDM, especially in models that allow for additional light relics (ACDM+Neg), based on 
the results for that model from [15]. 

Of course, the CMB power spectra also contain vastly more information than just 0x, 
and we also explore the consequences for cosmological models from the combination of this 
full CMB information with DESI BAO. We use as our baseline the temperature (TT) and po- 
larisation (EE) auto-spectra, plus their cross-spectra (TE), as incorporated in the Commander 
(for multipoles £ < 30) and plik (for £ > 30) likelihoods from the official PR3 release [160]. 
Subsequently, a new Planck data release PR4 has been made available, which involves a 
consistent reprocessing of the data from both the LFI and HFI instruments on Planck us- 
ing the new common pipeline NPIPE, leading to slightly more data, lower noise, and better 
consistency between frequency channels [161]. Various teams have released additional high-£ 


6This covariance is available from https://github.com/schoeneberg/2024 bbn. results. 
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likelihood packages using PRA [162-164], which are updated versions of likelihoods included 
in earlier Planck data releases, and use alternative methods to derive information from the 
high- power. However, in most scenarios we consider in this paper, results from these up- 
dated likelihoods are very similar to those obtained from the official PR3 plik likelihood. 
Therefore, we choose to keep plik as our baseline, and note any variations in the results due 
to differences between this and the CamSpec likelihood of [162, 163] in Appendix B." 

In addition to the primary temperature and polarization anisotropy power spectra, 
modern CMB experiments are also able to measure the power spectrum of the gravitational 
lensing potential, Coe, from the connected 4-point function [167, 168]. The latest and most 
precise CMB lensing data comes from the combination of NPIPE PR4 Planck CMB lensing 
reconstruction [169] and the Data Release 6 of the Atacama Cosmology Telescope (ACT) 
[22, 170, 171].8 We adopt this combined CMB lensing likelihood from both experiments as 
our baseline. For the sake of brevity, in the following text and figures, we will denote results 
obtained using temperature and polarisation information from Planck, and CMB lensing 
information from the Planck+ACT combination, simply as “CMB”. Where occasionally 
necessary, we will explicitly label results that do not use CMB lensing as “CMB (no lensing)”. 


2.4.3 Type Ia supernovae 


Type Ia supernovae (SN Ia) serve as standardizable candles which offer an alternative way 
to measure the expansion history of the universe. Historically, SN Ia led to the discovery of 
the accelerating expansion [1, 2], following earlier, more complex arguments for A-dominated 
models based on observations of large-scale structure [172, 173]. Within the ACDM model 
SN Ia have lower statistical power than modern BAO measurements, but provide important 
information on dark energy when analyzing the less restricted models considered in Section 5. 

In this paper, we make use of three different SN Ia datasets. The Pantheon+ compilation 
[174] consists of 1550 spectroscopically-confirmed SN Ia in the redshift range 0.001 < z < 
2.26. We use the public likelihood from [24] incorporating the full statistical and systematic 
covariance, imposing a bound z > 0.01 in order to mitigate the impact of peculiar velocities in 
the Hubble diagram [175]. More recently, the Union3 compilation of 2087 SN Ia, many (1363 
SN Ia) in common with Pantheon+, was presented in [25], and includes a different treatment 
of systematic errors and uncertainties based on Bayesian Hierarchical Modelling. Finally, the 
Dark Energy Survey, as part of their Year 5 data release, recently published results based on a 
new, homogeneously selected sample of 1635 photometrically-classified SN Ia with redshifts 
0.1 < z « 1.3, which is complemented by 194 low-redshift SN Ia (which are in common 
with the Pantheon+ sample) spanning 0.025 < z < 0.1 [26]. The contribution of SN Ia 
in general, and the differences between the different datasets!! are discussed in Section 5. 
In the rest of the paper, we will denote the Pantheon+ dataset as PantheonPlus, and the 


"CamSpec is bundled as an additional likelihood in the cobaya sampling code [165, 166], at https: //cobaya. 
readthedocs.io/en/latest/likelihood planck.html. 

*'The likelihood is available from https://github.com/ACTCollaboration/act, dr6. lenslike; we use the 
actplanck baseline option. 

?Data provided by the Union3 team, private communication. 

Data available at https://github.com/des-science/DES-SN5YR. 

11We note that these SN la data sets are not independent of each other and we therefore do not combine 
their results: Pantheon+ and Union3 share ~ 1360 supernovae but differ in their analysis methodology and 
marginalization over astrophysical and systematic parameters, while DES contributes a new high-z data set 
of ~ 1500 photometrically-classified SN la but still uses about 194 historical SN la at z < 0.1 in common with 
the other two. 
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parametrization parameter default prior 
background-only Qm — u (0.01, 0.99] 
no rq calibration rah (Mpc) — A10, 1000] 
with rq calibration Hp (kms^! Mpc!) — u [20, 100] 

Wp = U[0.005, 0.1 
CMB Wedm — u [0.001, 0.99] 

Wh — U[0.005, 0.1 

1000uc — û [0.5, 10] 

In(1019 A,) — u[1.61,3.91 

ns = U[0.8, 1.2] 

T — U (0.01, 0.8] 
extended OK 0 U[-0.3, 0.3] 

wo or w —1 U[—3, 1] 

Wa 0 U[—3, 2] 

>) my (eV) 0.06 U(0, 5] 

Neg 3.044  U[0.05,10] 


Table 2. Parameters and priors used in the analysis. All of the priors are flat in the ranges given. 
We consider two parametrizations, “background-only” when using BAO data only, and “CMB” where 
data from Planck and ACT are involved. In both cases, the same priors are used for parameter 
extensions. A single massive neutrino of mass )m, = 0.06 eV is assumed, except in the ACDM+m, 
model, for which we consider three degenerate massive neutrino species (N, = 3 in CAMB). In addition 
to the flat priors on wo and wa listed in the table, we also impose the requirement wo + Wa < 0 in 
order to enforce a period of high-redshift matter domination. 


DES-SN5YR dataset as DESY5, in order to provide a consistent, concise set of labels suitable 
for tables and figure legends, and to avoid ambiguities with the ‘+’ symbol used to denote 
the combinations of datasets. 


2.5 Cosmological inference 


We included in the cosmological inference code cobaya [165, 166] the PantheonPlus, Union3 
and DESY5 SN la likelihoods, as well as our new DESI BAO likelihoods. The CMB likeli- 
hoods used are based on public packages that are either included in the public cobaya version 
or available directly from the respective teams. 

When running cobaya, we rely on the Boltzmann code CAMB [176, 177] for theoretical 
cosmology calculations. When using the combined Planck+ACT lensing likelihood, we use 
higher precision settings as recommended by ACT. 

All Bayesian inference is performed using the Metropolis-Hastings MCMC sampler [178, 
179] in cobaya. Table 2 summarises the cosmological parameters that are sampled over in 
different runs, and the priors that are placed on them. In the base ACDM model, for data 
combinations that only probe the background evolution (i.e., BAO and SN Ia), we either 
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sample in the parameters Qm and rgh or—when using external data such as from BBN in 
order to help calibrate the rg value and break the rgh degeneracy in BAO data—in Om, Ho 
and the physical baryon density wp. When also including CMB likelihoods, we sample instead 
in the standard six-parameter basis (Weams Wb, 100010, In (1019 As), ns, T), where vc is an 
approximation to the acoustic angular scale 6,, A, is the amplitude of the primordial scalar 
power spectrum and n, is its spectral index, and 7 is the reionization optical depth. We 
also consider extensions of the minimal ACDM model in which other parameters are allowed 
to vary: the spatial curvature Ck, the sum of neutrino masses 5) mp, the single constant 
dark energy equation of state parameter w or the two parameters (wo, Wa) when allowing a 
time-varying equation of state w(a) [180, 181], and the effective number of relativistic species 
Neg. Priors on these extended parameters are shown in the bottom section of Table 2; in 
any instance where one of these parameters is not varied, it is held fixed at the default value 
listed. 

When including CMB data, we take advantage of the hierarchy between (fast) nuisance 
parameters and (slow) cosmological parameters by taking more steps along the latter, setting 
the oversampling parameter oversample power to 0.4. We also use the so-called dragging 
method that enables taking larger steps in the slow parameter subspace with the help of 
intermediate transitions of fast parameters [182]. For each dataset and model combination, 
we run four chains in parallel, starting from proposal covariance matrices built from pre- 
liminary runs. Chains are stopped when the Gelman-Rubin [183] criterion R — 1 < 0.01 is 
satisfied, where R is the largest eigenvalue of the ratio of the inter- to intra-chain covariance 
matrices. We further require the effective sample size of the chains to be > 10% to achieve 
percent precision on the moments of the marginal posteriors. In the case of symmetric 1D 
marginalized posteriors we report the mean and standard deviation of the samples, while 
the 68% minimal credible interval is quoted in other cases, except when otherwise stated. 
We use getdist!? [184] to derive the constraints presented in this paper. Where appropri- 
ate, we also compute the best-fit with the iminuit [185, 186] algorithm, as implemented in 
cobaya, starting from the four maximum a posteriori (MAP) points of the corresponding 
four MCMC chains. In the following, we note Aa — --2A log £ the difference (times —2) 
in log-posteriors at the maximum posterior points. 


3 DESI Distance Measurements 


3.1 Distance-Redshift Results 


Our summary Table 1 provides the BAO distance scales measured from each of the DESI 
tracers and redshift bins. As noted above, for the BGS and QSO samples, we report a 
single result for the angle-averaged quantity Dy/ra [79]. For all other redshift bins, we 
quote marginalized constraints on both Dy /ra and Dy/ra individually; note however that 
these two quantities are correlated with each other, with a correlation factor that depends 
on redshift, also provided in Table 1 and fully accounted for in the BAO likelihoods. 

Figure 1 shows a summary of these DESI BAO results in the form of a Hubble diagram. 
In order to conveniently display the results from all redshift bins together, we consistently 
convert from the (Dw/ra, Dy /rq) basis to (Dv /rq, Dy / Dg) for all tracers as labelled.!? The 
top row then shows Dy /rq (scaled by an arbitrary factor of 272/3 for visualisation purposes) 


https: //github.com/cmbant/getdist 
13In the (Dv /ra, Dwı / Du) basis the measurements from discrete tracers are much more independent of each 
other (though not totally so). The Lya BAO measurements still retain a significant correlation in this basis. 
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Figure 1. Top row: DESI measurements of the BAO distance scales at different redshifts, 
parametrized as (left) the ratio of the angle-averaged distance Dy = (zD?2,Dg)'/? to the sound 
horizon at the baryon drag epoch, ra, and (right) the ratio of transverse and line-of-sight comoving 
distances Fap = Dw/Dn, from all tracers and redshift bins as labeled. For visual clarity and to 
compress the dynamic range of the plot, an arbitrary scaling of 272/% has been applied on the left, 
and z-! on the right. The solid and dashed grey lines show model predictions from, respectively, the 
flat ACDM model that best fits this data, and from a ACDM model with parameters matching the 
Planck best-fit cosmology. The BGS and QSO data points appear only in the left panel and not the 
right one because the signal-to-noise ratio of the data is not yet sufficient to measure both parameters 
for these tracers. Bottom row: The same data points and models as in the top row, but now shown 
as the ratio relative to the predictions for the best-fit flat ACDM model. 


in the left panel, and Dy/Dg (similarly arbitrarily scaled by z71!) in the right panel. The 
solid and dashed grey lines in each panel indicate the corresponding model predictions for 
the ACDM model that best fit the DESI data (Section 4.1), and the Planck best-fit ACDM 
model, respectively. The lower panel shows the same data again but now as the ratio of the 
Dy/rq and Fap = Du/Dy values to those for the best-fit ACDM model to DESI data. The 
solid and dashed grey lines in these panels therefore represent the same two models as in the 
top row. 


3.2 Internal consistency of DESI results 


Figure 1 shows visually that the flat ACDM model provides a good fit to the DESI BAO 
results: quantitatively, the x? value for this fit is 12.66 for 10 degrees of freedom (dof), 
as we have 12 data points and 2 free parameters, namely Qm and Hora (Table 2). These 
two parameters have a direct relationship to the BAO data points shown in Figure 1, since 
in the flat ACDM model Nm fully determines FAp(z) and fixes the shape of Dy/ra as a 
function of redshift, while Horą sets a redshift-independent constant normalization term for 
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Figure 2. Left panel: 68% and 95% credible-interval contours for parameters Qm and rah obtained 
for a flat ACDM model from fits to BAO measurements from each DESI tracer type individually, as 
labeled. Results from all tracers are consistent with each other and the change in the degeneracy 
directions arises from the different effective redshifts of the samples. Right panel: the corresponding 
results in flat ACDM for fits to BAO results from all DESI redshift bins (blue), the final SDSS 
results from [139] (orange), and the combination of these two as described in the text (green). The 
corresponding result from the CMB (including CMB lensing) is shown in pink. 


Dy/ra. The single most anomalous result is the measurement of FAp(z = 0.51), which is 
mildly offset (at the ~ 2g significance level) from the ACDM expectation. This appears to 
be a simple statistical fluctuation and we do not regard it as significant. In Appendix A we 
explicitly confirm that this fluctuation does not alter any of the conclusions drawn about any 
cosmological models in Sections 4 to 7 below. 

It is also worth confirming that the parameters inferred from each individual redshift 
bin are consistent with each other. This is shown in Figure 2, where in the left panel 
we show the 68% and 95% credible intervals on Qm and Hora obtained from fitting a flat 
ACDM model to each DESI tracer type individually, using the priors described in Table 2. 
The slightly anomalous FAp(z — 0.51) value results in a mild shift in the contour for the 
0.4 < z < 0.6 LRG1 bin, but we do not regard this as significant. The individual contours 
have characteristic degeneracy directions that depend on redshift [187], and show a strong 
degree of overlap, confirming the internal consistency of the DESI data. A consequence of 
the change in the degeneracy directions with redshift is that the combination of all tracers 
provides a tight final constraint, shown in the right panel of Figure 2. 


3.3 Comparison to BAO results from SDSS 


'The region of the sky and the range of redshifts observed by DESI partially overlap with 
those observed by the previous generation of BOSS [188] and eBOSS [189] survey programs 
of SDSS [190], whose final BAO results were presented in [139]. It is therefore pertinent to 
compare our BAO results to those from SDSS**. As a rough guide, ~ 70% of the DESI DR1 
footprint was covered by BOSS; conversely, ~ 65% of the BOSS footprint has been covered 
by DESI DR1. Therefore the input catalog data used in the DESI and SDSS BAO analyses, 


14179] also shows a comparison to BAO measurements from 6dFGS [191] and WiggleZ [192], although these 
have lower precision 
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while different due to the details of the instrument performance and observing strategy, 
nevertheless have a fraction of shared objects, introducing a correlation that depends on the 
redshift and tracer type. 

This correlation is greatest at low redshifts, where the DESI BGS and LRG samples 
overlap the SDSS Main Galaxy Sample and LRGs from BOSS (LOWZ and CMASS) and 
eBOSS. There is also a substantial overlap in volume between the QSO samples of DESI and 
eBOSS in the redshift range 0.8 < z < 2.1. Despite this, however, fewer than 15% of the 
DESI quasars are also present in the eBOSS catalog. As the measurements of QSO power 
spectra are shot-noise dominated in both DESI and eBOSS, this results in a negligible overall 
correlation between the results. 

For redshifts below z = 0.6, the DESI DR1 data currently covers a smaller effective 
volume than SDSS (compare the V.g values in Table 1 to those reported in [139]), but this 
will change in future data releases. The redshift binning used differs between DESI and 
SDSS, which makes a direct comparison of BAO results complicated at z « 0.4, but the 
0.4 < z < 0.6 and 0.6 < z < 0.8 redshift bins closely match the redshift ranges used by BOSS 
and eBOSS.!? The SDSS results in these bins are 


Dyw/ra(z = 0.51) = 13.36 4 0.21, | SDSS LRG ai 

Dg/ra(z = 0.51) = 22.33 + 0.58, | 0.4 < z < 0.6, i 
and 

Du/ra(z = 0.70) = 17.86 + 0.33, | SDSS LRG T 

Dg/ra(z = 0.70) = 19.33 + 0.53, | 0.6 < z < 0.85, : 


respectively, which can be compared to the DESI results in Table 1. While the results at 
effective redshift z = 0.51 are in good agreement, a larger difference can be seen in the 
0.6 < z < 0.8 redshift bin, particularly in comparison to the DESI result Du /ra(z = 0.71) = 
16.85 + 0.32. 

To gauge the significance of this disagreement, we can estimate the degree of corre- 
lation between power spectrum measurements from DESI and SDSS from Eq. (2.7), using 
approximate values for the total volume overlap, the number density of galaxies in com- 
mon between the two surveys, and assuming a representative amplitude of the power spec- 
trum Po = 104 (h-  Mpc)? for both samples. This gives indicative values of C = 0.35 for 
0.4 < z < 0.6 and C = 0.21 for 0.6 < z < 0.8. Assuming that this degree of correlation 
also applies to BAO results derived from the power spectra, and accounting for the predicted 
changes in Dy /rq and Dg /rq between z = 0.70 and z = 0.71 in the DESI fiducial cosmology, 
the discrepancy between the DESI and SDSS results in the 0.6 « z « 0.8 redshift range is at 
roughly the ~ 30 level. 


The cause of this difference is not clear. As discussed in [79], numerous improvements 
in the data processing and BAO fitting methods have been introduced for DESI, and repro- 
cessing the raw BOSS/eBOSS catalog data using the DESI BAO pipeline gives small shifts 
compared to the results published in [139]. However, this accounts for a shift of at most 
a small fraction of the published uncertainties in Dy /rq and Dy/ra. We note that in this 


15eBOSS actually used the redshift range 0.6 < z < 1.0 to define the target sample, but the number density 
of LRGs in the catalog fell of very sharply beyond z > 0.85 so this approximation holds. 
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redshift bin the SDSS catalog was a composite formed from the BOSS CMASS galaxy sam- 
ple extending to z ^ 0.75 and the deeper eBOSS LRG sample extending to z = 1 over a 
much smaller sky footprint; in contrast the DESI LRG sample is much more uniform. The 
cause of the difference may simply an unlucky sample variance fluctuation—if so, this will 
soon become clear with later DESI data releases, which will have a much larger Veg in this 
redshift range. 

While the degree of correlation between BAO results from the galaxy and quasar samples 
for DESI and SDSS has only been approximately estimated as described above, in the Lya 
forest BAO analysis, the degree of correlation has been more thoroughly quantified in [80] 
and shown to agree very well. [80] also provides results for a combined “DESI+SDSS” Lya 
BAO measurement: 


Dy/ra(z = 2.33) = 38.80 + 0.75, (3.3) 
and 

Dg/ra(z = 2.33) = 8.72 + 0.14, (3.4) 
with an anticorrelation of p = —0.48 between Dw/rq and Dy/ra. 


The Dm/ra(z = 0.71) measurement from LRGs yields the greatest difference between 
SDSS and DESI at any redshift. Nevertheless, using the BAO measurements from all red- 
shifts together, the BAO distance ladders obtained from DESI and SDSS result in consistent 
inference of cosmological parameters. The right panel of Figure 2 shows the posterior cred- 
ible intervals for the parameters Qm and rq in a flat ACDM cosmology from fitting to the 
DESI (blue) and SDSS (orange) BAO data. The DESI values are described in Section 4 and 
Table 3. The results are clearly in very good agreement with each other, with no significant 
difference in Qm and a shift of just ~ lo in rah. 

In this paper our primary focus will be on the cosmological consequences of the DESI 
DR1 BAO results alone, as a homogeneously analysed dataset across redshift and tracer 
type. However, it is also possible to construct a dataset of BAO distance measurements 
from the combination of DESI and SDSS results in order to maximise measurement precision 
across the entire redshift range 0.1 < z < 4.2. Bearing in mind that the degree of correlation 
between BAO results from discrete galaxy and quasar tracers in the two surveys has not been 
precisely quantified—to avoid double-counting information—this combined sample should be 
selected by choosing the result from the survey covering the larger effective volume Veg at a 
given redshift. Thus, the composite BAO dataset can be constructed as follows: 


e at z < 0.6 where SDSS currently has a larger Ver, we use the SDSS results at zef = 
0.15,0.38 and 0.51 in place of the DESI BGS and lowest-redshift LRG points; 


e at z > 0.6 where DESI has Veg larger than that of SDSS, we use the DESI results from 
LRGs over 0.6 < z < 0.8, the LRG+ELG combination over 0.8 < z < 1.1, and ELGs 
and QSOs at higher redshifts; and 


e for the Lya BAO we use the combined DESI+SDSS result from Eqs. (3.3) and (3.4) 
above. 


We use the label “DESI+SDSS” to refer to this composite BAO dataset, while reiterating 
that this is not the same as simply combining the likelihoods from each survey individu- 
ally due to the overlap in volumes. The posterior in parameters (Qm, rah) inferred from 
fitting a flat ACDM model to this dataset is shown in green in the right panel of Figure 2. 
In Appendix A we present constraints using the DESI--SDSS data combination in various 
cosmological models, specifically comparing them to constraints using DESI data alone. 
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4 Cosmological constraints in the ACDM model 


4.1 Flat background 


In this section, we present cosmological constraints for the flat ACDM cosmological model in 
which Ok = 0, QA = 1— Qm. As discussed in Sections 3.2 and 3.3 above, in this model BAO 
distance scales alone constrain only two free parameters: the matter density parameter Qm 
and the combination Hora of the Hubble constant Hy and the sound horizon rg. The 68% 
credible-interval results for these parameters are 


Qm = 0.295 + 0.015, 
rah = (101.8 + 1.3) Mpc, 


) DESI BAO, (4.1) 


where to simplify units and notation we quote results for rah = Hora/(100 kms”! Mpe”!) 
instead of Hora. The posterior constraints are shown in the right panel of Figure 2, and it 
is clear that they are in very good agreement with the previously reported SDSS values of 
Qm = 0.299 + 0.016 and rah = (100.4 + 1.3) Mpc [139]. 

It is apparent from Figure 2 that DESI BAO prefer a somewhat larger value of rgh 
compared to the value rgh = (98.82 + 0.82) Mpc obtained from the combination of CMB 
temperature, polarisation and lensing (or rgh = (98.9 + 1.0) Mpc from the CMB when ex- 
cluding CMB lensing). It is therefore of interest to quantify the level of consistency between 
these two data sets. We do so by calculating the relative x? between the two datasets as 


x? = (pa — pp)! (Cova + Covg) + (pa — pg), (4.2) 


and converting it to the probability-to-exceed between two datasets (e.g. [193, 194]). Here 
Pa and pg refer to (Qm, rah) parameter vectors obtained from fits to the two datasets, and 
Cova and Covg to the corresponding 2 x 2 covariances. Applying this statistic to DESI 
DR1 BAO and CMB, we find that the differences in this parameter space amount to a mild 
1.9c-level discrepancy when including CMB lensing, or 1.60 without. We thus consider that 
DESI BAO data are consistent with the CMB in this parameter space. 

Since BAO distance measurements alone are sensitive to the combination Hora [195], 
an external calibration of the sound horizon rg is required in order to break the Ho-rq 
degeneracy and obtain a constraint on the Hubble constant Ho. This method of calibrating 
the BAO distance scales using the sound horizon at early times is known as the "inverse 
distance ladder" approach [196, 197]. Directly calibrating the BAO standard ruler using the 
value rq — 147.09 + 0.26 Mpc obtained from using all CMB and CMB lensing information 
[15] gives 


Ho = (69.29 + 0.87) kms! Mpe! (DESI BAO + rg from CMB). (4.3) 


An alternative approach that is more independent of CMB information is to use a prior on 
Qph? obtained from BBN (Section 2.4.1), which is sufficient to determine rg and break the 
degeneracy [68, 198]. Within flat ACDM, DESI DR1 BAO + BBN give 


Ho = (68.53 + 0.80) kms ! Mpe! (DESI BAO + BBN). (4.4) 


This constraint is 2096 more precise than the equivalent one obtained from SDSS BAO 
data, Hp = (67.35 + 0.97) kms”! Mpc^! [139], which is a result of the improved precision of 


—— = 


Ho 


model/dataset Ou 1090 % Or Wo Wa 
[kms¬t Mpc7!}] 
Flat ACDM 
DESI 0.295 + 0.015 = = = — 
DESL-BBN 0.295 + 0.015 68.53 + 0.80 — — — 
DESL-BBN--6, 0.2948 + 0.0074 68.52 + 0.62 — — — 
DESI+CMB 0.3069 + 0.0050 ^ 67.97 + 0.38 — — - 
ACDM+9k 
DESI 0.284 + 0.020 — 65163 — — 
DESI+BBN+9, 0.296 + 0.014 68.52 + 0.69 0.3168 — — 
DESI+CMB 0.3049 + 0.0051 68.51 Æ 0.52 2.4+1.6 — — 
wCDM 
DESI 0.293 + 0.015 BÊ = —0.99*0-15 — 
DESI--BBN-4-6, 0.295 + 0.014 08.0771 — —1.002* 0-021 — 
DESI+CMB 0.281 + 0.013 11312 = ES uror — 
DESI+CMB+Panth. 0.3095 ++ 0.0069 67.74 + 0.71 0.997 + 0.025 — 
DESI+CMB+Union3 0.3095 + 0.0083 67.76 + 0.90 0.997 + 0.032 -— 
DESI+CMB+DESY5 0.3169 +0.0065 66.92 + 0.64 0.967 + 0.024 — 
wogwg; CDM 
DESI 0.344+9:937 — — —0.55*039 < -1.32 
DESI+BBN+9, o. 65.0133 — —0.53*022 « —1.08 
DESI+CMB 0:344 7025 64.7135 — —0.45*031  —1.79*9108 
DESI+CMB+Panth. 0.3085 + 0.0068 68.03 + 0.72 0.827 + 0.063  —0.75*022 
DESI+CMB+Union3 0.3230 + 0.0095 66.53 + 0.94 0.65 +0.10 1.274538 
DESI+CMB+DESY5 0.3160 + 0.0065 67.24 + 0.66 0.727 + 0.067  —1.05*031 
wow; CDM+9k 
DESI 0.313 + 0.049 = Shae —0.70*0:22 < -1.21 
DESI+BBN+9, 03402005. 65.8139 5.0 08 0.5277 < -1.44 
DESI+CMB 0.34742: 931 64.32 -0.9+2 —-041*038 < -1.61 
DESI+CMB+Panth. 0.3084 +0.0067  68.06+0.74  0.3+18  —0.831+0.066  —0.73+9:%2 
DESL-CMB--Union3 0.3233*00059 66.45 + 0.98 0.4+1.9 0.64 + 0.11  —1.30*0$5 
DESI+CMB+DESY5 0.3163 ++ 0.0065 67.19 + 0.69 0.2 + 1.0 —0.725+0.071 —1.06*035 


Table 3. Cosmological parameter results from DESI DR1 BAO data in combination with external 
datasets and priors, in the baseline flat ACDM model and extensions including spatial curvature 
and two parametrizations of the dark energy equation of state, as listed. Results are quoted for the 
marginalized means and 68% credible intervals in each case, including for upper limits. Note that 
DESI data alone measures rgh and not Ho, but for reasons of space this result is omitted from the 
table and provided in the text instead. In this and other tables, the shorthand notation “CMB” is 
used to denote the addition of temperature and polarisation data from Planck and CMB lensing data 
from the combination of Planck and ACT. 
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Figure 3. Left panel: marginalized posterior constraints on matter density Qm and the Hubble 
constant Ho, obtained from combining DESI BAO data with external data used to calibrate the 
sound horizon rg, in a flat ACDM cosmological model. The combinations shown use a prior on wp 
determined from BBN (blue), the combination of a BBN wp, prior and measurement of the acoustic 
angular scale 9, (orange), and rg directly calibrated from CMB results from Planck (green). The 
pink contour shows the corresponding constraints from the combination of CMB and CMB lensing. 
Right panel: The marginalized 1D posteriors on Qm in flat ACDM, from DESI BAO, CMB and the 
three SN datasets, as labelled. 


the DESI BAO measurements at higher redshifts, leading to narrower posterior constraints 
perpendicular to the main degeneracy direction in the Qm-rah plane (Figure 2). The central 
value is also shifted by about ~ lo to higher Ho, a direct consequence of the similar shift in 
Hora shown in Figure 2. However, the two results are fully consistent with each other. 

We can also use a conservative model-independent prior on the acoustic angular scale 0, 
seen from the CMB (Section 2.4.2), without including additional CMB information. While 
on its own knowledge of 6, is insufficient to break the ra-^ degeneracy (because the baryon 
density Q5? remains unknown, so rq cannot be fixed), it is extremely robust and, in combi- 
nation with BBN information, further tightens the Ho result to 


Ho = (68.52 4 0.62) kms ! Mpc ! (DESI BAO+ BBN + 6.4), (4.5) 


a 0.9% measurement of the Hubble constant. As discussed further in Section 6 below, this 
more precise constraint can also be extended to be robust to many early-universe extensions of 
the base ACDM model, in particular to assumptions about the effective number of relativistic 
species Neg. 

It is interesting to note that the direct rg calibration from the CMB produces the 
weakest constraint on Ho despite apparently using the most external information and being 
the most model-dependent. This is because BBN and 0, information affect Qm and are 
applied consistently in the posterior sampling, and thus leverage the degeneracy between Qm 
and ah from BAO data, while a prior directly applied on rg cannot do so. 

The left panel of Figure 3 summarises the constraints in the Q.m-Ho plane obtained 
from the combination of DESI BAO data with each of the external priors discussed so far, 
and compares them to the combined CMB result from Planck and ACT. All combinations 
including DESI data favour somewhat higher values of Ho and lower values of Qm than the 
mean values for the CMB posterior. However, the results are not in serious tension: using 
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Figure 4. 68% and 95% marginalized posterior constraints on Qy—Q, plane (left) and Qm-Nx (right) 
in the one-parameter extension of the ACDM model with free curvature, ACDM+Qx. In the left panel 
the supernova contours are truncated at the lower-left by the U[—0.3, 0.3] prior on Ox. 


the metric outlined in Eq. (4.2) above, the biggest discrepancy is at the ~ 2.10 level, for 
DESI+BBN+9, compared to CMB. Given this level of agreement, there is no issue with 
combining DESI and CMB data to obtain joint constraints; doing so we find 


DESI BAO+ (4.6) 


Qm = 0.3069 + 0.0050, 
CMB. 


Ho = (67.97 + 0.38) kms”? Mpc”! 


These results are summarised in Table 3, which also shows parameter constraints obtained 
in other extended models. 

A final instructive comparison within the context of the flat ACDM is between the 
constraints on the matter density Nm offered by DESI and SN la. These are shown in the right 
panel of Figure 3. PantheonPlus reports Qm = 0.3314 0.018, Union3 gives Qm = 03595005. 
and DESY5 gives Qm = 0.353+0.017. There is therefore a moderate variation in both central 
values and quoted uncertainties across different SN Ia compilations, but all of them prefer 
higher values of Qm than DESI and the CMB. The statistical significance of the differences 
compared to DESI, calculated as described above, stands at 1.60 for PantheonPlus, 2.0c for 
Union3, and 2.60 for DESY5. While they do not meet a 3c threshold for significant tension, 
these numbers indicate a degree of disagreement between these datasets and DESI results 
when interpreted in the flat ACDM model. Should these mismatches persist and become 
more significant when more data is acquired, they will require further investigation. 


4.2 ACDM model with free spatial curvature 


Relaxing the condition of spatial flatness in the ACDM model, we allow the curvature pa- 
rameter Qx to vary. In an FLRW background, this is equivalent to allowing the dark energy 
density QA = 1 — O4, — Ox to vary independently from the matter density Qm, while still 
keeping the dark-energy equation of state fixed at w — —1. Because DESI provides rela- 
tive measures of the BAO scale at multiple redshifts, it can determine the expansion rate 
as a function of redshift and thus measure QA independent of any calibration of the sound 
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horizon from external data [199]. In this model DESI alone thus measures the parameters 
(Qm, Q4, rah), and finds the following 68% credible-interval constraints 


Om = 0.284 + 0.020, 
QA 005170055, DESI BAO. (4.7) 
rah = (100.9 + 1.6) Mpc, 


Expressed in terms of curvature, the DESI result is Ok = (0050095. Figure 4 shows the 


corresponding 2D credible intervals in both (Qm, Q4) and (Q4, Qg) plane. 

The addition of model-independent information on the acoustic angular scale 0, ef- 
fectively extends the redshift range of the BAO distance ladder, reaching all the way to 
last scattering at z ~ 1090. This dramatically improves the curvature constraints, giving 
OK = DOLOR Ee Adding a BBN prior this is further tightened to 


Ox = 0.0003170-0685 (DESI BAO + BBN + 6,). (4.8) 
The corresponding curvature constraint from the CMB is 


Qk = —0.0102 + 0.0054 (CMB). (4.9) 


This value is ~ 20 away from Qx = 0, and the difference is even larger (over 30) when not in- 
cluding CMB lensing. This is a well-known issue in the Planck PR3 likelihood (see discussion 
in [15], as well as, e.g., [200, 201]) caused by the combination of the geometric degeneracy 
limiting CMB measurements of curvature, and the observed slight excess of smoothing of 
high-/ peaks in the temperature power spectrum. Later Planck PRA likelihoods somewhat 
alleviate, but do not completely remove, this CMB preference for Oy < 0 [162-164]. Thus 
BAO results are crucial in providing a curvature constraint that is independent of the CMB. 
'The combined result from DESI BAO and all CMB information is 


Qk = 0.0024 + 0.0016 (DESI BAO+CMB). (4.10) 


5 Dark energy 


Over the past quarter century, an impressive variety of cosmological observations have con- 
firmed and vastly strengthened the evidence that the expansion of the universe is accelerating. 
The standard cosmological model has a flat background geometry, but a sub-critical total 
(cold and baryonic) matter density. While the measurements thus far are in good agreement 
with the simplest flat ACDM scenario—a vacuum energy described by the cosmological con- 
stant dominating the energy density and responsible for the late-time acceleration—much 
effort is dedicated to measuring the expansion and growth history of the universe and look- 
ing for any departures from this model. In this regard, constraining dark energy models 
beyond the simplest ACDM model is a principal goal of DESI. 
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Figure 5. Constraints on Qm and w in the flat CDM model. The constraints from DESI BAO 
alone are shown in blue, those from the CMB in pink, and different SN Ia compilations in solid and 
dashed green. The orange contour shows the combined constraint from DESI, CMB and PantheonPlus 
SN Ta. All contours show 68% and 95% credible intervals. Note the remarkable complementarity of 
cosmological probes in this plane. 


5.1 Flat wCDM model 


Although a cosmological constant fits existing data well, the tiny observed value of A relative 
to typical scales in particle physics poses great theoretical challenges [202-204]. Acceleration 
physics beyond A necessarily has dynamics — time dependence (and spatial perturbations, 
though these have diminishing effect the closer the time dependence is to A, i.e. constant). 

At the background cosmology level entering cosmic distances, the acceleration physics 
can be treated as an effective dark energy density and pressure. Thus we have a dark energy 
equation of state parameter, or pressure to energy density ratio, w(a), and a current dark 
energy density value, (qe, to describe the dark energy component. In the early 2000s, with the 
time variation w(a) inaccessible to observations, analyses often fixed w = const (wCDM), but 
this is insensitive to crucial dynamics that might be indicated by data. wCDM can however 
still be useful as an alert if the recovered constraint on constant w has statistically significant 
deviation from w = —1. Note the converse is not true: measuring w = —1 assuming w is 
constant does not indicate A is correct (known as the “mirage of A” [205]). 

Figure 5 shows the constraints on Qm and w from a variety of different data and com- 
binations. In this plane the DESI contours are close to vertically aligned, providing a tight 
constraint on Qm that is largely independent of w, in contrast to SN Ia and CMB probes, 
which show distinctive degeneracy directions corresponding to the transverse comoving dis- 
tances that each of these probes constrains [206]. We find 


Qm = 0.293 + 0.015, ê 
w = 0.994215, DESI BAO, 
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Figure 6. Left panel: 68% and 95% marginalized posterior constraints in the wy-w, plane for the 
flat wow; CDM model, from DESI BAO alone (black dashed), DESI + CMB (pink), and DESI + 
SN Ia, for the PantheonPlus [24], Union3 [25] and DESY5 [26] SNIa datasets in blue, orange and 
green respectively. Each of these combinations favours wy > —1, w, < 0, with several of them 
exhibiting mild discrepancies with ACDM at the = 20 level. However, the full constraining power is 
not realised without combining all three probes. Right panel: the 68% and 95% marginalized posterior 
constraints from DESI BAO combined with CMB and each of the PantheonPlus, Union3 and DESY5 
SN la datasets. The significance of the tension with ACDM (wo = —1, wa = 0) estimated from the 
Ax2,4p values is 2.50, 3.50 and 3.90 for these three cases respectively. 


from DESI alone, while combining DESI BAO with BBN and 0, significantly tightens the 


constraint on w to w = —1.002* 0.021. Adding CMB data shifts the contours slightly along 


the CMB degeneracy direction, giving 


Om = 0.281 + 0.013, 


= +0.062 
w = —1.122+0-062 


) DESI BAO+CMB. (5.2) 


Finally, the tightest constraints are obtained from the combination of these data with SN la. 
For example for the PantheonPlus SN la dataset: 


Qm = 0.3095 + 0.0069, ) DESI+CMB (5.3) 


w = —0.997 + 0.025, +PantheonPlus. 


Similar constraints are obtained when substituting PantheonPlus SN Ia for DESY5 or Union3 
(though with slightly larger uncertainties in the latter case). These results are summarised 
in Table 3. In summary, DESI data, both alone and in combination with other cosmological 
probes, do not show any evidence for a constant equation of state parameter different from 
—1 when a flat wCDM model is assumed. 


5.2 Flat wow; CDM model 


Taking into account the physical dynamics of dark energy, the parametrization w(a) = wo + 
Wa (1 — a) was derived and has been demonstrated to match the background evolution of 
distances arising from exact dark energy equations of motion to an accuracy of ~ 0.1% 


7 


for viable cosmologies over a wide range of physics — scalar fields, modified gravity, phase 
transitions [181, 207]. In this section, we present constraints on this model, referred to 
as wugwg; CDM, which reduces to ACDM for wo = —1, wa = 0. Constraining the y—Wa 
parameter space and its corresponding behavior as well as distinguishing it from A is a key 
science goal of DESI. 

We adopt wide flat priors on wo and wa (Table 2), together with the condition wo+Wa < 
0 imposed to enforce a period of high-redshift matter domination. Since the parameter space 
we explore includes models whose equation of state crosses the w — —1 boundary, we use the 
parametrized post-Friedmann approach [208] to compute the dark energy perturbations when 
calculating the CMB angular power spectrum. Figure 6 shows the marginalized posteriors in 
the wy-wa plane from DESI and combinations with other external datasets. DESI alone does 
not have sufficient power to break the degeneracy between wo and wa and thus the results 
are cut off by our priors (see Table 2), 


= +0.39 
wo = —0.55*039. 


5.4 
ND ) DESI BAO, (5.4) 


with the upper bound on wa referring to the 68% limit. This represents a mild pull away from 
the ACDM value, with a Axa rap between the maximum a posteriori of the wowaCDM model 
and the maximum of the posterior fixing (wo, wa) = (—1,0) of just —3.7 for 2 additional 
degrees of freedom. The cause of this preference for wọ > —1 is primarily due to the Fap 
measurement from the 0.4 < z < 0.6 LRG bin, which lies slightly higher (at the ~ 20 level) 
than the best-fit ACDM model can accommodate, as shown in Figure 1. In order to better fit 
this data point, the equation of state w(z S 0.5) of the best-fit wow; CDM model prefers to 
be high, thus preferentially pulling wo to less negative values than —1. On the other hand, to 
fit the other DESI points which are all fairly close to the ACDM. predictions, the parameter 
Wa prefers to be strongly negative in order to compensate the integrated effect of w(z) for 
those quantities at higher redshift. 

Given the small Ax? Ap it is clear there is no statistical preference for wow; CDM from 
DESI alone. CMB data alone also gives Ax?;4p = —3.7 for the MAP wowaCDM model 
compared to fixing (wo, wa) = (—1,0), again showing no statistical preference. Nevertheless, 
given the overlap of the CMB and DESI contours in the (wo > —1, wa < 0) quadrant, the 
combined results give 


wo = —0.45*034. 
Toas | DESIBAO+CMB (5.5) 
Wwa = —1.79*948, 


and the Axe, Ap decreases to —9.5, indicating a preference for an evolving dark energy equa- 
tion of state at the ~ 2.60 level. The contours in this scenario are however still impacted by 
the priors we have assumed, thus care is required in interpreting these shifts. 

SN Ia data alone allow for wa < 0 and, as shown in the left panel of Figure 6, in combina- 
tion with DESI BAO they also marginally favor wo > —1 although the statistical significance 
of this preference depends on the particular SN Ia dataset and is not overwhelming in any 
case. In order to break the degeneracy in the wo—w, plane it is necessary to look at the 
joint constraints from the combination of DESI, CMB and SN Ia probes, shown in the right 
panel of Figure 6. These constraints are now not prior-dominated in either parameter. We 
find that the results and the associated uncertainties again vary depending on the choice of 
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supernova dataset, particularly when comparing PantheonPlus to the other two. However, in 
all cases, the results are compatible with each other. We find marginalized posterior means 


Wa = -0.7516-26, +PantheonPlus, 
from combination with PantheonPlus, 
üa == 1277 bot. +Union3, 
with Union3 SN la, and 
Wa = -1.055031, +DESY5, 


when using the DESY5 SN Ia data. The Ax2;4p values between the maximum a posteriori of 
the woWaCDM model and the maximum of the posterior fixing (wo, Wa) = (—1,0) are —8.7, 
—15.2 and —18.1 for the combinations of DESI and CMB with PantheonPlus, Union3 and 
DESY5 respectively. These correspond to preferences for a wow; CDM model over ACDM 
model at the 2.50, 3.50 and 3.90 significance levels, respectively. ** 

A Bayesian model-selection analysis gives Bayes factors of | ln Baı| = 0.65, 2.4, and 
2.8 in favor of the wow; CDM model over ACDM for the combinations of DESI--CMB with 
PantheonPlus, Union3 and DESY5, respectively. On Jeffreys’ scale [209, 210], these indicate 
a weak preference for wow, CDM over ACDM by the first of these data combinations, and 
moderate preference by the latter two. 

For this analysis, we ran additional nested-sampling chains using the PolyChord sampler 
[211] from which we derived the Bayesian evidence and Bayes factors using the anesthetic 
software [212]. We have verified that the posteriors derived from these auxiliary chains 
match those derived from the main MCMC chains described in Section 2.5, and that the 
Bayes factors reported here provide consistent interpretations with those estimated from 
the main MCMC chains through the Savage-Dickey density ratio approximation [210]. Not 
unexpectedly given the current constraints, we find the exact values of the Bayesian evidence 
(hence the Bayes factors) mildly depend on the specific prior ranges for wo and wa (see 
Table 2). 

The Deviance Information Criterion (DIC) for model selection [213] indicates similar 
preferences: A(DIC) = DICwow,com — DICAcpM of —2.0, —7.4 and —7.7 for the same 
combinations, respectively (see Table 1 of [214] for a DIC reference scale). 

Further useful information in the wow, CDM analysis is given by the pivot redshift, 
Zp, and the pivot equation of state, wp = w(zp), which inform us about the redshift and 


16Based on several independent MAP estimates, the uncertainty on these significance levels is of order 0.20. 
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equation-of-state value at which w(z) is best constrained by a given data set or combination 
of data sets, see e.g. [215, 216]. We use here the formalism and definitions of [216] and find 


(0.57, 1.094 + 0.070)  (DESI+CMB) 


(0.26, —0.982 + 0.028)  (DESI-CMB+PantheonPlus) 
(zp, Wp) = (5.9) 
(0.33, —0.960 + 0.038) (DESI+CMB+WUnion3) 


(0.26, —0.946 + 0.026) (DESI+CMB+DESY5). 


Figure 7 provides a summary of the constraints on the expansion history from the 
flat cosmological models that we have explored. Here, we show the quantity ra/Du = 
H(z)ra as a function of redshift z, rescaled for convenience by 1 + z. The top, middle, 
and bottom panels respectively show derived constraints on the temporal evolution of the 
Hubble parameter in ACDM, wCDM and wow, CDM. models from fitting to DESI data, with 
the shaded regions indicating the 68% and 9596 credible regions. The same dashed line in 
all panels indicates the best fit Planck ACDM model. Because H(z)/(1 + z) is directly 
proportional to à, i.e. the expansion rate, its slope with redshift is directly proportional to 
minus the cosmic expansion acceleration, —à. The range of the variation of H(z) shows that 
in each of the three cosmological models, a period of accelerated expansion (ie. a negative 
slope) is required at z < 0.7. 


5.3 wgw;, CDM model with free spatial curvature 


Finally, we also explore the most general case of a time-varying equation of state w(a) together 
with allowing free spatial curvature, in the so-called *Og9w; CDM model”. Even in this 
broad parameter space, the combination of DESI, Planck and SN Ia is able to provide quite 
tight constraints. These are shown in Figure 8 for each of the PantheonPlus, Union3 and 
DESY5 SN Ia datasets. As seen before, the resultant central values and uncertainties vary 
depending on the choice of SN Ia data, but all three are still broadly consistent with each 
other, and all combinations place very tight constraints on the allowed range of deviations 
from spatial flatness. We find 


wo = —0.831 + 0.066, 
We = —0.73*032. DESHECMB (5.10) 


PantheonPI 
Oe 0:0003 +: 0.0018, J. en 
wo = —0.64 + 0.11, 
130023, DESI+CMB EN 
E "ٍ +Union3, 


Ok = —0.0004 + 0.0019, 


and 


wo = —0.725 + 0.071, 
iy == 106018), DESEHGMB (5.12) 


Oz = —0.0002 + 0.0019, | * PESYS. 


In comparison to the results for the flat wow; CDM model in Section 5.2, the addition of 
an extra degree of freedom in Qg leads to a slight broadening of the constraints in the 1y—wa 
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Figure 7. Model-dependent constraints on the redshift dependence of the Hubble parameter times 
the sound horizon, H(z)ra = ra/Du (scaled by 1/100/(1 + z) for visual clarity) in three different 
classes of models fit to all DESI data from Figure 1. From top to bottom, the panels show the 68% 
and 95% credible regions in the ACDM, wCDM and wgw, CDM. models respectively in the coloured 
bands. The dashed line in each panel shows the behaviour in the best fit Planck ACDM model. For 
convenience, the scale on the right-hand axis shows (ra/rbl?nck) H(z)/(1 + 2). 


plane, thus marginally reducing the significance of the tension with ACDM. Nevertheless, the 
trend observed remains the same, with all combinations of DESI + CMB + SN Ia preferring 
wo > —1 and wa < 0. 


6 Hubble constant 


The determination of the Hubble constant has been contentious for many decades [217— 
219]. By the turn of the century, a consensus value of around 70 kms”? Mpc^! had emerged 
[220]. However, since the first Planck results [221], a growing tension has emerged between 
Ho determinations based on physics of the early universe, which tend to cluster close to 
the Planck preferred value of 67 kms ! Mpc ^! (e.g., [15, 139]), and local distance-ladder 
measurements based on Cepheids or other anchors, which mostly prefer larger values around 
73 kms ! Mpc ^! (e.g., [222-225]). Although distance ladder measurements based on the 
Tip of the Red Giant Branch (TRGB) method [226] prefer a lower Ho than Cepheid-based 
calibrations, they currently have larger uncertainties, hence do not yet provide a conclusive 
assessment. This tension between the CMB and Cepheid-based determinations stands at 
the ~ 4-50 level in ACDM and, if not due to some unidentified residual systematics, may 
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Figure 8. Marginalized posteriors on wo, Wa and (kx in a model with a time-varying dark energy 
equation of state and free spatial curvature, from DESI and CMB data combined with SN Ia from 
PantheonPlus, Union3 and DESY5 in blue, orange and green respectively. All combinations provide 
tight limits on Ox. Constraints on wo and wa in each case broaden a little compared to those shown 
in the flat case (Figure 6) but the overall trend remains the same. 


indicate deficiencies in the standard cosmological model. For reviews of the Hp tension, see 
[227—229]. 

As discussed in Sections 2.4.1 and 4.1, BAO alone cannot provide an absolute distance 
measurement, but can constrain Ho when combined with some external information capable 
of calibrating the sound horizon scale and breaking the rg—h degeneracy. Within the flat 
ACDM cosmological model, calibration of the DESI DR1 BAO with primordial deuterium 
abundances and BBN and the acoustic angular scale measurement leads to very precise Ho 
determinations summarised in Eqs. (4.4) and (4.5), the tightest of which is Ho = 68.52 + 
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Figure 9. 68% credible-interval constraints on the Hubble constant, assuming the flat ACDM model. 
The blue, bold whiskers show DESI BAO measurements in combination with an external BBN prior 
on 9h? and measurement of the acoustic angular scale 0,, the BBN prior alone, or with the CMB 
measurement of the sound horizon, rT'^"**, The thin blue whiskers show the corresponding results 
from SDSS BAO and the combination of DESI+SDSS BAO results, as labelled. The orange whiskers 
show the results from CMB anisotropy measurements from Planck and ACT, while the green whiskers 
show measurements of Ho from the distance ladder with either Cepheids or TRGB. 
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Figure 10. Left panel: The 68% and 95% credible-interval contours for Qm and Hora obtained from 
fitting DESI DR1 BAO data in the base flat ACDM model and in four extension models which modify 
the background geometry or late-time expansion history. Right panel: A summary of the tension in 
the Ho measurements obtained from the DESI BAO results combined with other data, and the SHOES 
result of [225], assuming different cosmological models. 


L6 = 


0.62 kms”! Mpc”!. As shown in Figure 9, while this DESI4+BBN+9, constraint is slightly 
higher than the Planck CMB value, it remains in tension (at the 3.70 level) with the SHOES 
Cepheid-based distance-ladder result Ho = 73.04 + 1.04 kms”! Mpc^! (see Figure 10). 

When the assumption of the flat ACDM model is relaxed, the o constraints obtained 
when adding DESI BAO measurements (in combination with external priors) are loosened 
too. There are two qualitatively different ways in which the more general cosmological models 
allow more freedom and lead to weaker H constraints. 

First, the calibration of the sound horizon using BBN relies on assumptions about the 
physics at the time of BBN (as well as on the CMB temperature To, which however is 
measured very accurately [230]). In particular, if the effective number of relativistic degrees 
of freedom Neg is allowed to vary from its value Neg = 3.044 in the base ACDM model, the 
value of wp inferred from light element abundances is altered. When accounting for this in 
the wp prior as discussed in Section 2.4.1 using the results of [156], but otherwise keeping the 
late-time geometry unchanged, the Ho result from Eq. (4.5) relaxes slightly to 


Ho = (68.5 + 1.4) kms ! Mpc ! (DESI + BBN +6,4, free Neg), (6.1) 


a 2% measurement. Keeping the BBN prior but dropping the 6, information changes the 
direction of the contours in the 2D (Qm, Ho) parameter space but the marginalized 1D con- 
straint remains very similar, Ho = (68.5 + 1.5) kms”! Mpc^!. When adding the full CMB 
information to DESI BAO, this result slightly shifts to 


Ho = (68.3 + 1.1) kms ! Mpc} (DESI BAO+CMB, free Ng). (6.2) 


'The second type of additional freedom that leads to weakened H constraints is allowed 
by beyond-ACDM models which change the geometry or the late-time expansion of the 
universe. This directly affects the BAO measurement of rgh, thus degrading the precision 
with which Ho can be determined, even when the early-universe physics remains the same. 
The left panel of Figure 10 shows the joint constraint in the Qm- rah plane as the amount 
of freedom in the parameterization of the expansion history is increased (for the models 
considered in the previous section). Clearly, the uncertainty on rgh increases dramatically 
as more model freedom is allowed. Nevertheless, the central values do not move much, and 
the results remain consistent with each other across models: we find rgh = (101.9 € 1.3) Mpc 
(ACDM), rah = (101.0 + 1.6) Mpe (ACDM4-Ox), rah = (101.7132) Mpc (wCDM), rah = 
(96.4733) Mpe (wow; CDM) and rah = (98.2152) Mpc (wowaCDM+Qx). As a consequence, 
the Ho values determined from combining DESI BAO with BBN or BBN--0, priors also 
remain stable or slightly decrease, as shown in Table 3. This conclusion remains true in most 
cases when combining DESI BAO with the CMB, or with the CMB and SN la information 
(the one exception is in the (CDM model, where the CMB pulls to higher Ho in the absence 
of SN Ia), i.e., independent of the late-time modifications studied, BAO data in combination 
with other probes always prefer low central values of Ho. However, as the Hp uncertainties 
increase in these models, the level of tension with SHOES still decreases, as summarised in 
the right panel of Figure 10. Further results on the value of Hp in models which combine 
extensions to the background expansion history with changes to the neutrino sector are listed 
in Table 4. 


7 Neutrinos 


A generic prediction of the hot Big Bang model is a relic neutrino background which leaves 
detectable imprints on cosmological observations. Neutrinos are the only known particles to 
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Figure 11. Left panel: The marginalized 1D posterior constraints on X` m, from different combina- 
tions of datasets, in the single parameter extension flat ACDM 4-5 m, model. As explained in the 
text, here we use a model with 3 degenerate mass eigenstates and with a minimal prior $5 m, > 0 eV. 
The minimal masses for the normal or inverted mass ordering scenarios correspond to > m, > 0.059 
eV and Y m, > 0.10 eV respectively, shown by the vertical dashed lines and the shaded regions. 
Right panel: Joint marginalized 68% and 95% credible intervals on $` m, and Ho from Planck, CMB 
and DESI+CMB data, illustrating the degeneracy between these parameters from the CMB, and how 
DESI BAO data contribute to breaking it. The vertical shaded regions indicate the 68% and 95% con- 
straints on Ho from DESI BAO data combined with knowledge of 6, and Qh? in the ACDM+) m, 
model. This shows how DESI BAO breaks the primarily geometric degeneracy to place an upper limit 


on )) mp. 


behave as radiation in the early universe and as dark matter at late times, so they affect both 
the acoustic oscillations in the primordial plasma as well as the background evolution and 
structure formation. Cosmological observations are sensitive to both the number of neutrino 
species and their total mass (e.g., [69]), making cosmology constraints complementary to 
terrestrial neutrino experiments. 


7.1 Sum of neutrino masses 


The base model we have adopted so far assumes the sum of neutrino masses to be $5 m, = 
0.06 eV, with a single massive eigenstate and two massless ones. This is motivated by the 
lower bound for X` m, from neutrino oscillation experiments. In this section, we consider a 
single-parameter extension beyond this minimal model in which m, is allowed to freely 
vary, in order to explore the constraining power on ` m, of DESI data. Amongst terrestrial 
experiments directly measuring the neutrino mass, KATRIN [231] has produced the tightest 
constraints to date, from measuring the endpoint of the tritium P-decay spectrum. This 
gives an upper bound on the effective electron-neutrino mass that is independent of the 
cosmological model of mg < 0.8eV (90% CL) [232], equivalent to » 5 m, < 2.4eV (90% CL). 

Neutrino oscillation experiments have also shown that at least two of the three active 
neutrino masses are non-zero, but the ordering of these masses is not known. In the normal 
ordering or normal hierarchy (NH), the two lowest mass neutrino eigenstates have the smallest 
mass splitting, implying that the total neutrino mass must be )`) m, > 0.059 eV, while in the 
inverted hierarchy (IH), however, the smallest mass splitting occurs between the two highest 
mass eigenstates, necessitating a total mass of ` m, > 0.10eV [233]. When allowing X` m, 


— 35 = 


model / dataset Qu Ho [kms ! Mpc!] Xm, [eV] Nett 
ACDM+)Y m, 

DESI+CMB 0.3037 + 0.0053 68.27 + 0.42 < 0.072 — 
ACDM+Nog 

DESI+CMB 0.3058 + 0.0060 68.3 x 1.1 — 3.10 + 0.17 
wCDM-4- m, 

DESI+CMB 0.282 + 0.013 ram « 0.123 — 
DESI+CMB+Panth. 0.3081 + 0.0067 67.81 + 0.69 < 0.079 — 
DESI+CMB+Union3 0.3090 + 0.0082 67.72 + 0.88 < 0.078 — 
DESI4+CMB+DESY5 0.3152 + 0.0065 67.01 + 0.64 « 0.073 — 
wCDM+ Neg 

DESI+CMB 0.281 + 0.013 T1.011-8 — 2.97 + 0.18 
DESI+CMB+Panth. 0.3090 + 0.0068 67.9+1.1 — 3.07 +0.18 
DESI+CMB+Union3 0.3097 + 0.0084 67.8 1.2 — 3.06 + 0.18 
DESI+CMB+DESY5 0.3163 + 0.0067 67.2: 1.1 — 3.09 + 0.18 
WoWaCDM+)> m, 

DESI--CMB 0.344+0-032 64.7*21 « 0.195 
DESI+CMB+Panth. 0.3081 + 0.0069 68.07 + 0.72 < 0.155 — 
DESI+CMB+Union3 0.3240 + 0.0098 66.48 + 0.94 < 0.185 — 
DESI4+CMB+DESY5 0.3165 + 0.0069 67.22 + 0.66 « 0.177 — 
WoWa CDM +Nef 

DESI+CMB Dade. DS — 2.89 + 0.17 
DESI+CMB+Panth. — 0.3093 + 0.0069 67.5 X 1.1 — 2.93 + 0.18 
DESI+CMB+Union3 0.3245 + 0.0098 65.9 +1.3 — 2.91 + 0.18 
DESI+CMB+DESY5 0.3172 0.0067 66.6 + 1.1 — 2.92 + 0.18 


Table 4. Cosmological parameter results from DESI DR1 BAO data in combination with external 
datasets when considering extensions to the baseline ACDM model in the neutrino sector. Results 
with two-sided error bars refer to the marginalized means and 68% credible intervals; upper bounds on 
Y m, refer to the 95% limits. Note that the label “CMB” includes CMB lensing from the combined 
Planck+ACT likelihood. All constraints on )`) m, assume a model with 3 degenerate mass eigenstates 
and a minimal prior $, m, > 0 eV. (See Eq. (7.3) and Eq. (7.4) for the results using other priors.) 
The empty )) m, and Neg fields indicate that the fixed respective values of X` m, = 0.06eV and 


Neg = 3.044 


were adopted. 
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to be a free parameter, we adopt a model with three degenerate neutrino mass eigenstates. 
This degenerate mass model does not exactly correspond to either of the physically expected 
NH or IH scenarios; however, it produces a very good approximation of the observable effects 
of both [234]. In the event that a positive detection of non-zero neutrino mass is possible, 
an analysis using the degenerate mass model with a prior ))m,„ > 0 will recover the correct 
value of $5 m, for both NH and IH scenarios with little reconstruction bias [235]. On the 
other hand, when a positive detection is not possible, using the degenerate model with 
appropriately modified priors on ` m, as above will also recover the correct upper bounds 
for both the NH and IH scenarios [236]. 

'The sum of the neutrino masses affects cosmology in two ways. First, the high-velocity 
dispersion of neutrinos implies that they free-stream over large distances, thus suppressing 
the late-time clustering and power spectrum of matter at small scales, below the free stream- 
ing length. However, if X` m, is varied with other cosmological parameters along the CMB 
degeneracy, the net effect of neutrino masses is an overall almost-scale-independent suppres- 
sion of the amplitude of the matter power spectrum (together with a very small shift in the 
BAO scale). Second, at the low redshifts of relevance to DESI, neutrinos are non-relativistic 
(© m, > T ~ 107? eV) and therefore contribute to the total non-relativistic matter density 
Wm = Ut + We + wy, where Q, = )) m„/(93.14 eVh?), see e.g. [234]. Y m, thus affects the 
background evolution and in particular the redshift of matter-A equality. Forthcoming DESI 
analyses of the full shape of the galaxy power spectrum are sensitive to neutrino masses 
through the first effect, and constraints from them will be described in [81, 237]. The BAO 
data used in this paper only distinguish the background geometry. However, BAO constraints 
on the expansion history, when combined with the CMB and CMB lensing, are very helpful 
in improving neutrino mass constraints. This is because the CMB is sensitive to the neutrino 
mass via the angular diameter distance to recombination and via the effects of >) m, on 
the CMB lensing, and both of these effects can be mimicked in the CMB by varying other 
cosmological parameters such as Ho and wm. 

Thus while BAO data are not directly sensitive to the suppression effects of neutrinos 
on the power spectrum, by determining the background geometry at low redshifts they help 
to break the CMB degeneracy and constrain Ho, thus greatly tightening the upper limit 
on X` mp. This effect is illustrated in the right panel of Figure 11. In extensions of flat 
ACDM, the power of BAO to pin down the late-time expansion history diminishes, resulting 
in notably weaker >> m, constraints. 

For the Planck CMB alone, and assuming flat ACDM, the constraints are 


Som, <0.21eV (95%, CMB), (7.1) 


when including CMB lensing from the combination of Planck and ACT, while adding the 
DESI BAO data sharply reduces this to 


Y my <0.072eV (95%, DESI BAO+CMB). (7.2) 


The left panel of Figure 11 shows the corresponding 1D marginalized posteriors. The posterior 
peaks at }>m, = 0eV, which is excluded by terrestrial oscillation experiments, but the 
minimal mass ` m, = 0.059 eV is not excluded by the cosmological fits. We find a AX, AP 
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of 3.8 for the scenario with X` m, = 0.059eV compared to » m, = 0eV from DESI BAO 
and CMB. 

The result in Eq. (7.2) represents a substantial tightening of the upper bound compared 
to previous state-of-the-art results [22, 42, 238-241], of which the previous tightest reported 
upper bound was ` m, < 0.082 eV [240]. As can be seen from the right panel of Figure 11,the 
improvement relative to CMB-only constraints is driven primarily by the tighter constraints 
on Ho (and consequently Nm) obtained from BAO. The shift relative to the previous tightest 
upper bounds is instead a consequence of both the preference of DESI data for slightly higher 
rah and thus higher Ho (which is itself coupled to lower Qm), together with the improved 
lensing constraint from including the latest ACT (DR6) lensing data. Should the addition 
of more data in the future pull towards lower Ho, we may expect the neutrino mass limit 
to relax, even if the data achieve a higher precision. For the same reason, as discussed in 
Appendix A, the upper bound obtained using the combined BAO data from DESI and SDSS 
(Section 3.3) instead of DESI alone is also slightly looser, X` m, < 0.082eV, due to the 
slightly lower rgh and Hg values seen in Figures 2 and 9. 

We also caution that the upper bound obtained is strongly dependent on the choice of 
prior for 3; m,. While we have deliberately chosen to use Nm, > 0eV for our primary 
analysis, one can also impose physically motivated priors corresponding specifically to the 
NH Om, > 0.059 eV) and IH (>`) m, > 0.1eV) scenarios. Applying these alternative priors 
we find the following upper limit for NH: 


$m, > 0.059 eV), 


and for IH: 


Y^m,«0145ev (95%, DESI BAO+CMB; (7.4) 
>> my > 0.10 eV). 


Therefore cosmological results alone do not yet strongly favour the NH over the IH: the 
preference for NH over IH is at the 2g level. 

It is possible to combine the cosmology result for X` m, with constraints from oscillation 
and P-decay experiments to compute the Bayesian evidence in favour of one hierarchy over 
the other, although the results depend strongly on the choice of prior for the individual 
masses. A full calculation is beyond the scope of this paper, but even with a minimally 
informative objective Bayesian prior the cosmological constraints we find—assuming the flat 
ACDM background cosmology—when combined with terrestrial data give the Bayes factor 
for the normal hierarchy over the inverted one above 100 [242-244]. 

However, we repeat the caveat that the limits are substantially relaxed in more extended 
dark energy models that affect the background geometry, such as those in Section 5. Table 4 
shows the corresponding upper bounds for $` m, when allowing for a wCDM or wow; CDM 
background; as expected, these are significantly weaker: for instance, in a ugwa CDM. back- 
ground, the upper limit from DESI BAO and CMB relaxes to }>m, < 0.195 eV (for the 
>) my, > 0 eV prior). Importantly though, while allowing a time-varying DE equation of state 
parameter significantly loosens the upper bound on >> m,, allowing >) m, to vary freely only 
marginally affects the constraints on wo and Wa reported in Section 5.2 and therefore does 
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not substantially change the conclusions of that section. For example, when allowing ` m, 
to vary freely, we find 


= +0.068 
wo = —0.725 9076; DESI+CMB+DESY5, (7.5) 
Wa = —1.07*038, free ) mp 


which can be compared with Eq. (5.8). Combinations with the other two SN Ia datasets are 
similarly close to the results in Eqs. (5.6) and (5.7) obtained with fixed >`) m, = 0.06 eV. 


7.2 Number of effective relativistic species 


Under the standard assumption that three active neutrinos thermalize in the early universe, 
additional dark relativistic degrees of freedom can be parametrized in terms of a change to 
the effective number of neutrino species, Neg. This is defined such that the total relativistic 
energy density, after the annihilation of electrons with positrons, is given by 


4/3 

pv = Neg : (=) Pr» (7.6) 
where p~ is the photon energy density. In the standard cosmological model with three massive 
species of neutrinos and no other particles that are relativistic at recombination (other than 
photons), Neg — 3.044. Extended models with light sterile neutrinos or other dark relics 
(such as “dark radiation") that are generated well before recombination produce effects very 
similar to that of active neutrinos and so can be usefully explored in terms of constraints on 
Neg in a ACDM--Neg model. 

Similar to the case of >) m,, constraints on Neg from the CMB alone exhibit a geomet- 
rical degeneracy because changing the relativistic energy density before recombination shifts 
the sound horizon, and this effect can be absorbed through changes to Ho (or Qm), such 
that increasing Neg corresponds to higher Ho or lower Qm. Therefore the DESI BAO mea- 
surements at lower redshifts again contribute by breaking the geometric degeneracy through 
their ability to constrain Qm. 

In particular, while the combination of CMB anisotropies and lensing power spectra 
from Planck and ACT give 


Neg = 2.98 + 0.20 (CMB), (7.7) 


the addition of DESI BAO changes this to 


Neg - 3.10 + 0.17 (CMB+DESI BAO), (7.3) 


i.e., a small shift of the central value due to the DESI preference for lower Qm, and a ~ 15% 
reduction in the uncertainty. As shown in Table 4, despite the additional freedom allowed 
in the neutrino sector, the recovered value of Ho remains compatible with that from Planck. 
This remains true even when allowing the temporal variation of the dark energy equation of 
state. 

Note that changes in Neg can also produce damping of the BAO amplitude and shifts to 
the scale and phase of the BAO oscillations due to the neutrino dragging effect [245, 246]. In 
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principle these effects could introduce biases in the BAO measurements, if the template used 
for the fitting is drawn from a cosmological model that does not contain the same physical 
effects. The implications of this for BAO measurements are studied in [247], where it is 
found that even for large values Neg = 3.70 that are excluded at high significance by the 
joint constraints above, the systematic offset in the recovered BAO scale Dy /ra is, at most, 
~ 0.196. 


8 Conclusions 


In this paper we have presented the first cosmological results from DESI, the first Stage- 
IV galaxy survey in operation, marking the start of a new era of dark energy experiments. 
These results are based on samples of bright galaxies, LRGs, ELGs, quasars and Lya forest 
tracers in the redshift range 0.1 « z « 4.2. These data include a total of over 6 million 
unique extragalactic spectroscopic redshifts from the first year of DESI observations alone 
(out of the full five-year survey program), already representing an increase of more than a 
factor of two over the number of tracers used in the previous largest such dataset assembled 
[139], which was the culmination of two decades of observations with SDSS. This remarkable 
achievement is made possible by the speed and quality of the DESI instrument [70], whose 
increased spectral resolution also delivers significantly improved redshift accuracy. 

'The results in this paper are based on the measurement of the BAO scale using different 
DESI tracers of the matter density (galaxies, quasars, and the Lyman-a forest) in seven 
redshift bins, presented in detail in [79, 80]. The observed data were analysed using a state- 
of-the-art blind analysis pipeline to protect against confirmation bias and included end-to-end 
validation. This methodology enabled us to unleash the potential of the BAO method to, in 
conjunction with other probes, powerfully constrain cosmological parameters. 

Having demonstrated the internal consistency of the DESI BAO measurements over the 
full redshift range of observations and shown that the results are in agreement with previous 
measurements from SDSS, we proceeded to examine the constraints and implications of these 
data for a range of cosmological models. 

In the standard flat ACDM cosmology, we determine the matter density parameter 
Om = 0.295 + 0.015, and the product of the drag-epoch sound horizon and the scaled Hubble 
constant of rgh = (101.8 + 1.3) Mpc. These values are slightly different from those measured 
from the combination of CMB anisotropies from Planck plus CMB lensing from Planck and 
ACT, giving a somewhat lower value of Qm and a higher rgh, although the discrepancy is not 
statistically significant. In combination with a conservative prior on the baryon abundance 
from BBN, DESI BAO data determine the Hubble constant value to be Ho = (68.53 + 
0.80) kms”! Mpc”!, the most precise measurement to date that does not rely on information 
from CMB anisotropies. Combining the conservative BBN prior with an equally conservative 
prior on the extremely precise and model-independent measurement of the acoustic angular 
scale 0, gives Ho = (68.52 + 0.62) kms^! Mpc^!, approaching the precision from Planck. 
Each of these two H constraints is in a > 3g tension with the SHOES Cepheid-based distance- 
ladder result. In combination with the full CMB information, DESI results give Ho = 
(67.97 + 0.38) km s^! Mpc~!—a 0.6% precision measurement. 


BAO distance measurements are particularly important in constraining model exten- 
sions to ACDM, where they help break geometric degeneracies that limit the power of the 
CMB. We have examined several models of dark energy, allowing spatial curvature to vary. 
Together with CMB information, DESI BAO data provide extremely tight limits on the spa- 
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tial curvature in all scenarios, with Qx = —0.0024 + 0.0016 in the simplest single-parameter 
extension of ACDM and almost equally tight constraints in more extended model variations. 

'The measurement of the dark energy equation of state w is a key science goal of DESI. 
Assuming the wCDM model where the equation-of-state parameter is constant in time, we 
find w — —0.99*015 from DESI alone, and w = —0.997 + 0.025 from the combination of 
DESI BAO, CMB, and SN Ia results from the Pantheon+ compilation, in good consistency 
with ACDM. This result does not appreciably change when the Pantheon+ data are replaced 
by those from other recent SN Ia releases from Union3 and the Dark Energy Survey (DES- 
SN5YR). 

However, when the equation of state is allowed to vary with time, w(a) = wo+(1—a)wa, 
DESI data favour solutions with wg > —1 and wa < 0. The combination of DESI and CMB 
gives wo = —0.45*031 and Wa = —1.79*078, and indicates a ~ 2.20 difference to ACDM. 
When adding information from SN la, all combinations prefer wọ > —1 and wa < 0, with 
the level of the tension with ACDM remaining at the ~ 2.50 level for combining DESI and 
CMB information with Pantheon+, but increasing to 3.50 and 3.90 levels for the Union3 
and DESYSSN SN Ia datasets respectively. All three of the primary BAO, CMB and SN Ia 
probes contribute partially to this tension and the results including the 3 SN Ia datasets 
are mutually consistent with each other. Moreover, combining any two of the DESI BAO, 
CMB or SN data sets shows some level of departure from the ACDM model. Relaxing 
the assumption of a spatially flat geometry through varying Oy marginally increases the 
uncertainties but does not change the overall picture. It remains important to thoroughly 
examine unaccounted-for sources of systematic uncertainties or inconsistencies between the 
different datasets that might be contributing to these results. Nevertheless, these findings 
provide a tantalizing suggestion of deviations from the standard cosmological model that 
motivate continued study and highlight the potential of DESI and other Stage-IV surveys to 
pin down the nature of dark energy. 


Neutrinos are the only particles of the Standard Model of particle physics whose mass 
parameters are unknown. DESI, and Stage-IV surveys more generally, will improve cosmo- 
logical constraints on neutrino mass parameters and provide key insights into their mass 
hierarchy. Allowing the sum of the neutrino masses to vary, the combination of DESI BAO 
and CMB information breaks a geometric degeneracy (in CMB constraints) between the 
Hubble constant Ho and the amplitude of matter fluctuations and thus places an extremely 
tight upper bound on the total mass, $m, < 0.072eV (95% CL). We however caution 
that this substantial tightening of the upper bound compared to the previous state-of-the- 
art measurements is partly driven by the preference of DESI data for a higher value of Ho. 
Moreover, the upper bound obtained depends on the priors chosen for X` m, and, because 
DESI contributes through breaking a geometrical degeneracy, the upper limit is relaxed in 
extended models that alter the background geometry. We examined this point explicitly 
in the paper by considering changes to the upper bound on ` m, in models with a vary- 
ing dark energy equation of state. Finally, we have also reported an updated constraint 
Neg = 3.10 + 0.17 on the effective number of extra relativistic degrees of freedom from the 
combination of DESI BAO and CMB data, and showed that this constraint only marginally 
shifts to Neg = 2.89 + 0.17 even when varying the dark energy equation of state. 


As the first set of cosmological results from DESI, this paper, together with the ac- 
companying results in [79, 80], demonstrates the enormous power of the DESI instrument 
and survey. Subsequent papers will examine the implications of measurement of the full 
clustering broadband shape of DESI [81, 237] and constraints on primordial non-Gaussianity 
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[248]. With data quality and the speed of survey completion continuing to match or exceed 
expectations, future data releases will soon be able to provide even better insights into the 
hints of the exciting cosmological findings presented here. 


9 Data Availability 


The data used in this analysis will be made public along the Data Release 1 (details in 
https://data.desi.lbl.gov/doc/releases/). 


Acknowledgments 


This material is based upon work supported by the U.S. Department of Energy (DOE), Office 
of Science, Office of High-Energy Physics, under Contract No. DE-ACO02-05CH11231, and 
by the National Energy Research Scientific Computing Center, a DOE Office of Science User 
Facility under the same contract. Additional support for DESI was provided by the U.S. Na- 
tional Science Foundation (NSF), Division of Astronomical Sciences under Contract No. AST- 
0950945 to the NSF National Optical-Infrared Astronomy Research Laboratory; the Science 
and Technology Facilities Council of the United Kingdom; the Gordon and Betty Moore Foun- 
dation; the Heising-Simons Foundation; the French Alternative Energies and Atomic Energy 
Commission (CEA); the National Council of Humanities, Science and Technology of Mexico 
(CONAHCYT); the Ministry of Science and Innovation of Spain (MICINN), and by the DESI 
Member Institutions: https://www.desi.lbl.gov/collaborating-institutions. 

The DESI Legacy Imaging Surveys consist of three individual and complementary 
projects: the Dark Energy Camera Legacy Survey (DECaLS), the Beijing-Arizona Sky Sur- 
vey (BASS), and the Mayall z-band Legacy Survey (MzLS). DECaLS, BASS and MzLS 
together include data obtained, respectively, at the Blanco telescope, Cerro Tololo Inter- 
American Observatory, NSF NOIRLab; the Bok telescope, Steward Observatory, University 
of Arizona; and the Mayall telescope, Kitt Peak National Observatory, NOIRLab. NOIRLab 
is operated by the Association of Universities for Research in Astronomy (AURA) under a 
cooperative agreement with the National Science Foundation. Pipeline processing and analy- 
ses of the data were supported by NOIRLab and the Lawrence Berkeley National Laboratory. 
Legacy Surveys also uses data products from the Near-Earth Object Wide-field Infrared Sur- 
vey Explorer (NEOWISE), a project of the Jet Propulsion Laboratory /California Institute of 
Technology, funded by the National Aeronautics and Space Administration. Legacy Surveys 
was supported by: the Director, Office of Science, Office of High Energy Physics of the U.S. 
Department of Energy; the National Energy Research Scientific Computing Center, a DOE 
Office of Science User Facility; the U.S. National Science Foundation, Division of Astronom- 
ical Sciences; the National Astronomical Observatories of China, the Chinese Academy of 
Sciences and the Chinese National Natural Science Foundation. LBNL is managed by the 
Regents of the University of California under contract to the U.S. Department of Energy. 
The complete acknowledgments can be found at https: //www.legacysurvey.org/. 

Any opinions, findings, and conclusions or recommendations expressed in this material 
are those of the author(s) and do not necessarily reflect the views of the U.S. National Science 
Foundation, the U.S. Department of Energy, or any of the listed funding agencies. 

The authors are honored to be permitted to conduct scientific research on lolkam Du'ag 
(Kitt Peak), a mountain with particular significance to the Tohono O'odham Nation. 


- 42- 


References 


1 


10 


11 


12 


13 


14 


15 


16 


17 


SUPERNOVA SEARCH TEAM collaboration, Observational evidence from supernovae for an 
accelerating universe and a cosmological constant, Astron. J. 116 (1998) 1009 
[astro-ph/9805201]. 


SUPERNOVA COSMOLOGY PROJECT collaboration, Measurements of Q and A from 42 high 
redshift supernovae, Astrophys. J. 517 (1999) 565 [astro-ph/9812133]. 


J. Frieman, M. Turner and D. Huterer, Dark Energy and the Accelerating Universe, Ann. Rev. 
Astron. Astrophys. 46 (2008) 385 [0803.0982]. 


D.H. Weinberg, M.J. Mortonson, D.J. Eisenstein, C. Hirata, A.G. Riess and E. Rozo, 
Observational Probes of Cosmic Acceleration, Phys. Rept. 530 (2013) 87 [1201 . 2434]. 


SUPERNOVA SEARCH TEAM collaboration, Type la supernova discoveries at z > 1 from the 
Hubble Space Telescope: Evidence for past deceleration and constraints on dark energy 
evolution, Astrophys. J. 607 (2004) 665 [astro-ph/0402512]. 


ESSENCE collaboration, Observational Constraints on the Nature of the Dark Energy: First 
Cosmological Results from the ESSENCE Supernova Survey, Astrophys. J. 666 (2007) 694 
[astro-ph/0701041]. 


SNLS collaboration, Supernova Constraints and Systematic Uncertainties from the First 3 
Years of the Supernova Legacy Survey, Astrophys. J. Suppl. 192 (2011) 1 [1104.1443]. 


SUPERNOVA COSMOLOGY PROJECT collaboration, The Hubble Space Telescope Cluster 
Supernova Survey: V. Improving the Dark Energy Constraints Above z>1 and Building an 
Early- Type- Hosted Supernova Sample, Astrophys. J. T46 (2012) 85 [1105.3470]. 


M. Betoule, R. Kessler, J. Guy, J. Mosher, D. Hardin, R. Biswas et al., Improved cosmological 
constraints from a joint analysis of the SDSS-II and SNLS supernova samples, A&A 568 
(2014) A22 [1401.4064]. 


C.L. Bennett, A.J. Banday, K.M. Gorski, G. Hinshaw, P. Jackson, P. Keegstra et al., 
Four-Year COBE DMR Cosmic Microwave Background Observations: Maps and Basic 
Results, ApJ 464 (1996) L1 [astro-ph/9601067]. 


A.D. Miller, R. Caldwell, M.J. Devlin, W.B. Dorwart, T. Herbig, M.R. Nolta et al., A 
Measurement of the Angular Power Spectrum of the Cosmic Microwave Background from L — 
100 to 400, ApJ 524 (1999) L1 [astro-ph/9906421]. 


S. Hanany et al., MAXIMA-1: A Measurement of the cosmic microwave background 
anisotropy on angular scales of 10 arcminutes to 5 degrees, Astrophys. J. Lett. 545 (2000) L5 
[astro-ph/0005123]. 

BOOMERANG collaboration, A measurement by Boomerang of multiple peaks in the angular 
power spectrum of the cosmic microwave background, Astrophys. J. 571 (2002) 604 
[astro-ph/0104460). 

G. Hinshaw, D. Larson, E. Komatsu, D.N. Spergel, C.L. Bennett, J. Dunkley et al., Nine-year 


Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Parameter 
Results, ApJS 208 (2013) 19 [1212.5226]. 


Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi 
et al., Planck 2018 results. VI. Cosmological parameters, A&A 641 (2020) A6 [1807 .06209]. 


Z. Hou et al., Constraints on Cosmology from the Cosmic Microwave Background Power 
Spectrum of the 2500 deg? SPT-SZ Survey, Astrophys. J. 782 (2014) 74 [1212.6267]. 


S. Aiola, E. Calabrese, L. Maurin, S. Naess, B.L. Schmitt, M.H. Abitbol et al., The Atacama 
Cosmology Telescope: DR4 maps and cosmological parameters, J. Cosmology Astropart. Phys. 
2020 (2020) 047 [2007 .07288]. 


- 43 - 


[18] 


[19] 


[20 


[21] 


[22] 


[23] 


[24 


[25] 


[26 


[27] 


[28 


[29 


30 


31 


32 


33 


[34] 


C. Heymans, T. Tróster, M. Asgari, C. Blake, H. Hildebrandt, B. Joachimi et al., KiDS-1000 
Cosmology: Multi-probe weak gravitational lensing and spectroscopic galaxy clustering 
constraints, A&A 646 (2021) A140 [2007 . 15632]. 


T.M.C. Abbott, M. Aguena, A. Alarcon, S. Allam, O. Alves, A. Amon et al., Dark Energy 
Survey Year 3 results: Cosmological constraints from galaxy clustering and weak lensing, 
Phys. Rev. D 105 (2022) 023520 [2105 . 13549]. 


S. More, S. Sugiyama, H. Miyatake, M.M. Rau, M. Shirasaki, X. Li et al., Hyper 
Suprime-Cam Year 3 results: Measurements of clustering of SDSS-BOSS galaxies, 
galazy-galazy lensing, and cosmic shear, Phys. Rev. D 108 (2023) 123520 [2304. 00703]. 


Dark Energy Survey and Kilo-Degree Survey Collaboration, T.M.C. Abbott, M. Aguena, 
A. Alarcon, O. Alves, A. Amon et al., DES Y3 + KiDS-1000: Consistent cosmology 
combining cosmic shear surveys, The Open Journal of Astrophysics 6 (2023) 36 [2305 . 17173]. 


M.S. Madhavacheril, F.J. Qu, B.D. Sherwin, N. MacCrann, Y. Li, I. Abril-Cabezas et al., The 
Atacama Cosmology Telescope: DR6 Gravitational Lensing Map and Cosmological 
Parameters, ApJ 962 (2024) 113 [2304 . 05203]. 


Z. Pan, F. Bianchini, W.L.K. Wu, P.A.R. Ade, Z. Ahmed, E. Anderes et al., Measurement of 
gravitational lensing of the cosmic microwave background using SPT-3G 2018 data, 
Phys. Rev. D 108 (2023) 122005 [2308 . 11608]. 


D. Brout, D. Scolnic, B. Popovic, A.G. Riess, A. Carr, J. Zuntz et al., The Pantheon+ 
Analysis: Cosmological Constraints, ApJ 938 (2022) 110 [2202.04077]. 


D. Rubin, G. Aldering, M. Betoule, A. Fruchter, X. Huang, A.G. Kim et al., Union Through 
UNITY: Cosmology with 2,000 SNe Using a Unified Bayesian Framework, arXiv e-prints 
(2023) arXiv:2311.12098 [2311 .12098]. 


DES collaboration, The Dark Energy Survey: Cosmology Results With ^1500 New 
High-redshift Type Ia Supernovae Using The Full 5-year Dataset, ApJ (accepted, 2024) 
[2401 . 02929]. 


P.J.E. Peebles, Statistical Analysis of Catalogs of Extragalactic Objects. I. Theory, ApJ 185 
(1973) 413. 


P.J.E. Peebles and M.G. Hauser, Statistical Analysis of Catalogs of Extragalactic Objects. TII. 
The Shane- Wirtanen and Zwicky Catalogs, ApJS 28 (1974) 19. 


E.J. Groth and P.J.E. Peebles, Statistical analysis of catalogs of extragalactic objects. 7. Two 
and three point correlation functions for the high - resolution Shane- Wirtanen catalog of 
galaxies, Astrophys. J. 217 (1977) 385. 


M. Davis and P.J.E. Peebles, A Survey of galaxy redshifts. 5. The Two point position and 
velocity correlations, Astrophys. J. 267 (1982) 465. 


V. de Lapparent, M.J. Geller and J.P. Huchra, The Mean Density and Two-Point Correlation 
Function for the CfA Redshift Survey Slices, ApJ 332 (1988) 44. 


C.M. Baugh and G. Efstathiou, The three-dimensional power spectrum measured from the 
APM galaxy survey - I. Use of the angular correlation function., MNRAS 265 (1993) 145. 


J.A. Peacock, S. Cole, P. Norberg, C.M. Baugh, J. Bland-Hawthorn, T. Bridges et al., A 
measurement of the cosmological mass density from clustering in the 2dF Galaxy Redshift 
Survey, Nature 410 (2001) 169 [astro-ph/0103143]. 


W.J. Percival, C.M. Baugh, J. Bland-Hawthorn, T. Bridges, R. Cannon, S. Cole et al., The 
2dF Galaxy Redshift Survey: the power spectrum and the matter content of the Universe, 
MNRAS 327 (2001) 1297 [astro-ph/0105252]. 


— 44 — 


35 


36 


37 


38 


39 
40 


41 


42 


43 


44 


45 


46 


47 


48 


49 


50 


51 


[52] 


K. Abazajian, J.K. Adelman-McCarthy, M.A. Agüeros, S.S. Allam, S.F. Anderson, J. Annis 
et al., The First Data Release of the Sloan Digital Sky Survey, AJ 126 (2003) 2081 
[astro-ph/0305492]. 


M. Tegmark, M.R. Blanton, M.A. Strauss, F. Hoyle, D. Schlegel, R. Scoccimarro et al., The 
Three-Dimensional Power Spectrum of Galaxies from the Sloan Digital Sky Survey, ApJ 606 
(2004) 702 [astro-ph/0310725]. 


M. Tegmark, D.J. Eisenstein, M.A. Strauss, D.H. Weinberg, M.R. Blanton, J.A. Frieman 
et al., Cosmological constraints from the SDSS luminous red galaxies, Phys. Rev. D 74 (2006) 
123507 [astro-ph/0608632]. 


D.J. Eisenstein and W. Hu, Baryonic Features in the Matter Transfer Function, ApJ 496 
(1998) 605 [astro-ph/9709112]. 


D.J. Eisenstein, Dark energy and cosmic sound [review article], New A Rev. 49 (2005) 360. 


B. Bassett and R. Hlozek, Baryon acoustic oscillations, in Dark Energy: Observational and 
Theoretical Approaches, P. Ruiz-Lapuente, ed., p. 246 (2010), DOI. 


D.H. Weinberg, M.J. Mortonson, D.J. Eisenstein, C. Hirata, A.G. Riess and E. Rozo, 
Observational probes of cosmic acceleration, Phys. Rep. 530 (2013) 87 [1201 2434]. 


S. Alam, M. Aubert, S. Avila, C. Balland, J.E. Bautista, M.A. Bershady et al., Completed 
SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological implications from 
two decades of spectroscopic surveys at the Apache Point Observatory, Phys. Rev. D 103 
(2021) 083533 [2007 .08991]. 


C. Blake, D. Parkinson, B. Bassett, K. Glazebrook, M. Kunz and R.C. Nichol, Universal 
fitting formulae for baryon oscillation surveys, MNRAS 365 (2006) 255 [astro-ph/0510239]. 


H.-J. Seo and D.J. Eisenstein, Improved Forecasts for the Baryon Acoustic Oscillations and 
Cosmological Distance Scale, ApJ 665 (2007) 14 [astro-ph/0701079]. 


R.E. Angulo, C.M. Baugh, C.S. Frenk and C.G. Lacey, The detectability of baryonic acoustic 
oscillations in future galaxy surveys, MNRAS 383 (2008) 755 [astro-ph/0702543]. 


P. McDonald and D.J. Eisenstein, Dark energy and curvature from a future baryonic acoustic 
oscillation survey using the Lyman-a forest, Phys. Rev. D 76 (2007) 063009 
[astro-ph/0607122). 


M. McQuinn and M. White, On estimating Lya forest correlations between multiple sightlines, 
MNRAS 415 (2011) 2257 [1102.1752]. 


A. Font-Ribera, D. Kirkby, N. Busca, J. Miralda-Escudé, N.P. Ross, A. Slosar et al., 
Quasar-Lyman a forest cross-correlation from BOSS DR11: Baryon Acoustic Oscillations, J. 
Cosmology Astropart. Phys. 2014 (2014) 027 [1311.1767]. 


D.J. Eisenstein, I. Zehavi, D.W. Hogg, R. Scoccimarro, M.R. Blanton, R.C. Nichol et al., 
Detection of the Baryon Acoustic Peak in the Large-Scale Correlation Function of SDSS 
Luminous Red Galaxies, ApJ 633 (2005) 560 [astro-ph/0501171]. 


S. Cole, W.J. Percival, J.A. Peacock, P. Norberg, C.M. Baugh, C.S. Frenk et al., The 2dF 
Galaxy Redshift Survey: power-spectrum analysis of the final data set and cosmological 
implications, MNRAS 362 (2005) 505 [astro-ph/0501174]. 


W.J. Percival, S. Cole, D.J. Eisenstein, R.C. Nichol, J.A. Peacock, A.C. Pope et al., 
Measuring the Baryon Acoustic Oscillation scale using the Sloan Digital Sky Survey and 2dF 
Galaxy Redshift Survey, MNRAS 381 (2007) 1053 [0705.3323]. 


W.J. Percival, B.A. Reid, D.J. Eisenstein, N.A. Bahcall, T. Budavari, J.A. Frieman et al., 
Baryon acoustic oscillations in the Sloan Digital Sky Survey Data Release 7 galaxy sample, 
MNRAS 401 (2010) 2148 [0907 . 1660]. 


— 45 — 


[53] 


[54 


[55] 


[56 


[57] 


[58] 


[59] 


[60 


[61] 


[62 


[63] 


[64] 


[65] 


[66] 


[67] 


C. Blake, T. Davis, G.B. Poole, D. Parkinson, S. Brough, M. Colless et al., The WiggleZ Dark 
Energy Survey: testing the cosmological model with baryon acoustic oscillations at z— 0.6, 
MNRAS 415 (2011) 2892 [1105.2862]. 


C. Blake, E.A. Kazin, F. Beutler, T.M. Davis, D. Parkinson, S. Brough et al., The WiggleZ 
Dark Energy Survey: mapping the distance-redshift relation with baryon acoustic oscillations, 
MNRAS 418 (2011) 1707 [1108.2635]. 


E.A. Kazin, J. Koda, C. Blake, N. Padmanabhan, S. Brough, M. Colless et al., The WiggleZ 
Dark Energy Survey: improved distance measurements to z — 1 with reconstruction of the 
baryonic acoustic feature, MNRAS 441 (2014) 3524 [1401 .0358]. 


F. Beutler, C. Blake, M. Colless, D.H. Jones, L. Staveley-Smith, L. Campbell et al., The 6dF 
Galaxy Survey: baryon acoustic oscillations and the local Hubble constant, MNRAS 416 
(2011) 3017 [1106.3366]. 


P. Carter, F. Beutler, W.J. Percival, C. Blake, J. Koda and A.J. Ross, Low redshift baryon 
acoustic oscillation measurement from the reconstructed 6-degree field galaxy survey, MNRAS 
481 (2018) 2371 [1803 .01746]. 


L. Anderson, E. Aubourg, S. Bailey, D. Bizyaev, M. Blanton, A.S. Bolton et al., The 
clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: baryon 
acoustic oscillations in the Data Release 9 spectroscopic galaxy sample, MNRAS 427 (2012) 
3435 [1203.6594]. 


L. Anderson, É. Aubourg, S. Bailey, F. Beutler, V. Bhardwaj, M. Blanton et al., The 
clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: baryon 
acoustic oscillations in the Data Releases 10 and 11 Galaxy samples, MNRAS 441 (2014) 24 
[1312.4877]. 


S. Alam, M. Ata, S. Bailey, F. Beutler, D. Bizyaev, J.A. Blazek et al., The clustering of 
galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological 
analysis of the DR12 galaxy sample, MNRAS 470 (2017) 2617 [1607 . 03155]. 


M. Ata, F. Baumgarten, J. Bautista, F. Beutler, D. Bizyaev, M.R. Blanton et al., The 
clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 quasar 
sample: first measurement of baryon acoustic oscillations between redshift 0.8 and 2.2, 
MNRAS 473 (2018) 4773 [1705 .06373]. 


J.E. Bautista, R. Paviot, M. Vargas Magana, S. de la Torre, S. Fromenteau, H. Gil-Marín 

et al., The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: 
measurement of the BAO and growth rate of structure of the luminous red galaxy sample from 
the anisotropic correlation function between redshifts 0.6 and 1, MNRAS 500 (2021) 736 
[2007 . 08993]. 


J. Hou, A.G. Sánchez, A.J. Ross, A. Smith, R. Neveux, J. Bautista et al., The completed 
SDSS-IV extended Baryon Oscillation Spectroscopic Survey: BAO and RSD measurements 
from anisotropic clustering analysis of the quasar sample in configuration space between 
redshift 0.8 and 2.2, MNRAS 500 (2021) 1201 [2007 . 08998]. 


N.G. Busca, T. Delubac, J. Rich, S. Bailey, A. Font-Ribera, D. Kirkby et al., Baryon acoustic 
oscillations in the Lya forest of BOSS quasars, A&A 552 (2013) A96 [1211.2616]. 


T. Delubac, J.E. Bautista, N.G. Busca, J. Rich, D. Kirkby, S. Bailey et al., Baryon acoustic 
oscillations in the Lya forest of BOSS DR11 quasars, A&A 574 (2015) A59 [1404.1801]. 


J.E. Bautista, N.G. Busca, J. Guy, J. Rich, M. Blomqvist, H. du Mas des Bourboux et al., 
Measurement of baryon acoustic oscillation correlations at z = 2.3 with SDSS DR12 
Lya-Forests, A&A 603 (2017) A12 [1702.00176]. 


H. du Mas des Bourboux, J. Rich, A. Font-Ribera, V. de Sainte Agathe, J. Farr, 


— 46 — 


68 


69 


70 


71 


72 


73 


74 


75 


76 


TT 


78 


79 


80 


81 


82 


83 


84 


85 


86 


T. Etourneau et al., The Completed SDSS-IV Extended Baryon Oscillation Spectroscopic 
Survey: Baryon Acoustic Oscillations with Lya Forests, ApJ 901 (2020) 153 [2007 . 08995]. 


É. Aubourg, S. Bailey, J.E. Bautista, F. Beutler, V. Bhardwaj, D. Bizyaev et al., 
Cosmological implications of baryon acoustic oscillation measurements, Phys. Rev. D 92 
(2015) 123516 [1411.1074]. 


DESI Collaboration, A. Aghamousa, J. Aguilar, S. Ahlen, S. Alam, L.E. Allen et al., The 
DESI Experiment Part I: Science, Targeting, and Survey Design, arXiv e-prints (2016) 
arXiv:1611.00036 [1611.00036]. 


DESI Collaboration, B. Abareshi, J. Aguilar, S. Ahlen, S. Alam, D.M. Alexander et al., 
Overview of the Instrumentation for the Dark Energy Spectroscopic Instrument, AJ 164 
(2022) 207 [2205 . 10939]. 


A.P. Cooper, S.E. Koposov, C. Allende Prieto, C.J. Manser, N. Kizhuprakkat, A.D. Myers 
et al., Overview of the DESI Milky Way Survey, ApJ 947 (2023) 37 [2208 08514]. 


J. Moon, D. Valcin, M. Rashkovetskyi, C. Saulder, J.N. Aguilar, S. Ahlen et al., First 
detection of the BAO signal from early DESI data, MNRAS 525 (2023) 5406 [2304 . 08427]. 


DESI Collaboration, A.G. Adame, J. Aguilar, S. Ahlen, S. Alam, G. Aldering et al., The 
Early Data Release of the Dark Energy Spectroscopic Instrument, arXiv e-prints (2023) 
arXiv:2306.06308 [2306 . 06308]. 


K. Koyama, Cosmological Tests of Modified Gravity, Rept. Prog. Phys. 79 (2016) 046902 
[1504 . 04623]. 


A. Joyce, L. Lombriser and F. Schmidt, Dark Energy Versus Modified Gravity, Ann. Rev. 
Nucl. Part. Sci. 66 (2016) 95 [1601.06133]. 


M. Ishak, Testing General Relativity in Cosmology, Living Rev. Rel. 22 (2019) 1 [1806 . 10122]. 


S. Alam, C. Arnold, A. Aviles, R. Bean, Y.-C. Cai, M. Cautun et al., Towards testing the 
theory of gravity with DESI: summary statistics, model predictions and future simulation 
requirements, J. Cosmology Astropart. Phys. 2021 (2021) 050 [2011 05771]. 


DESI Collaboration, DESI 2024 II: Sample definitions, characteristics and two-point 
clustering statistics, in preparation (2024) . 


DESI Collaboration, A.G. Adame, J. Aguilar, S. Ahlen, S. Alam, D.M. Alexander et al., 
DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars, arXiv e-prints 
(2024) arXiv:2404.03000 [2404 . 03000]. 


DESI Collaboration, A.G. Adame, J. Aguilar, S. Ahlen, S. Alam, D.M. Alexander et al., 
DESI 2024 IV: Baryon Acoustic Oscillations from the Lyman Alpha Forest, arXiv e-prints 
(2024) arXiv:2404.03001 [2404 . 03001]. 


DESI Collaboration, DESI 2024 V: Analysis of the full shape of two-point clustering statistics 
from galaxies and quasars, in preparation (2024) . 


DESI Collaboration, DESI 2024 I: Data Release 1 of the Dark Energy Spectroscopic 
Instrument, in preparation (2025) . 


S. Brieden, H. Gil-Marín and L. Verde, A tale of two (or more) h's, J. Cosmology Astropart. 
Phys. 2023 (2023) 023 [2212.04529]. 


C. Blake and K. Glazebrook, Probing Dark Energy Using Baryonic Oscillations in the Galaxy 
Power Spectrum as a Cosmological Ruler, ApJ 594 (2003) 665 [astro-ph/0301632]. 


H.-J. Seo and D.J. Eisenstein, Probing Dark Energy with Baryonic Acoustic Oscillations from 
Future Large Galaxy Redshift Surveys, ApJ 598 (2003) 720 [astro-ph/0307460]. 


E.V. Linder, Baryon oscillations as a cosmological probe, Phys. Rev. D 68 (2003) 083504 
[astro-ph/0304001]. 


— 47- 


87 


88 
89 


90 


91 


92 


93 


94 


95 


96 


97 


102 


103 


104 


K. Thepsuriya and A. Lewis, Accuracy of cosmological parameters using the baryon acoustic 
scale, J. Cosmology Astropart. Phys. 2015 (2015) 034 [1409.5066]. 


D.J. Schlegel et al., lsdr9, in prep. (2023) . 


B. Flaugher, H.T. Diehl, K. Honscheid, T.M.C. Abbott, O. Alvarez, R. Angstadt et al., The 
Dark Energy Camera, AJ 150 (2015) 150 [1504.02900]. 


Dark Energy Survey Collaboration, T. Abbott, F.B. Abdalla, J. Aleksié, S. Allam, A. Amara 
et al, The Dark Energy Survey: more than dark energy - an overview, MNRAS 460 (2016) 
1270 [1601 . 00329]. 


H. Zou, X. Zhou, X. Fan, T. Zhang, Z. Zhou, J. Nie et al., Project Overview of the 
Beijing-Arizona Sky Survey, PASP 129 (2017) 064101 [1702.03653]. 


A. Dey, D. Rabinowitz, A. Karcher, C. Bebek, C. Baltay, D. Sprayberry et al., Mosaic3: a 
red-sensitive upgrade for the prime focus camera at the Mayall 4m telescope, in Ground-based 
and Airborne Instrumentation for Astronomy VI, C.J. Evans, L. Simard and H. Takami, eds., 
vol. 9908 of Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 
p. 99082C, Aug., 2016, DOI. 


C. Hahn, M.J. Wilson, O. Ruiz-Macias, S. Cole, D.H. Weinberg, J. Moustakas et al., The 
DESI Bright Galaxy Survey: Final Target Selection, Design, and Validation, AJ 165 (2023) 
253 [220808512]. 


R. Zhou, B. Dey, J.A. Newman, D.J. Eisenstein, K. Dawson, S. Bailey et al., Target Selection 
and Validation of DESI Luminous Red Galaxies, AJ 165 (2023) 58 [2208.08515]. 


A. Raichoor, J. Moustakas, J.A. Newman, T. Karim, S. Ahlen, S. Alam et al., Target 
Selection and Validation of DESI Emission Line Galaxies, AJ 165 (2023) 126 [2208.08513]. 


E. Chaussidon, C. Yéche, N. Palanque-Delabrouille, D.M. Alexander, J. Yang, S. Ahlen et al., 
Target Selection and Validation of DESI Quasars, ApJ 944 (2023) 107 [220808511]. 


A.D. Myers, J. Moustakas, S. Bailey, B.A. Weaver, A.P. Cooper, J.E. Forero-Romero et al., 
The Target-selection Pipeline for the Dark Energy Spectroscopic Instrument, AJ 165 (2023) 
50 [2208 . 08518]. 


J.H. Silber, P. Fagrelius, K. Fanning, M. Schubnell, J.N. Aguilar, S. Ahlen et al., The Robotic 
Multiobject Focal Plane System of the Dark Energy Spectroscopic Instrument (DESI), AJ 165 
(2023) 9 [2205 09014]. 


Raichoor et al., in preparation (2024) . 


DESI Collaboration, A.G. Adame, J. Aguilar, S. Ahlen, S. Alam, G. Aldering et al., 
Validation of the Scientific Program for the Dark Energy Spectroscopic Instrument, AJ 167 
(2024) 62 [2306 .06307]. 


T.-W. Lan, R. Tojeiro, E. Armengaud, J.X. Prochaska, T.M. Davis, D.M. Alexander et al., 
The DESI Survey Validation: Results from Visual Inspection of Bright Galaxies, Luminous 
Red Galaxies, and Emission-line Galaxies, ApJ 943 (2023) 68 [2208 08516]. 


D.M. Alexander, T.M. Davis, E. Chaussidon, V.A. Fawcett, A. X. Gonzalez-Morales, 
T.-W. Lan et al., The DESI Survey Validation: Results from Visual Inspection of the Quasar 
Survey Spectra, AJ 165 (2023) 124 [2208.08517]. 


E.F. Schlafly, D. Kirkby, D.J. Schlegel, A.D. Myers, A. Raichoor, K. Dawson et al., Survey 
Operations for the Dark Energy Spectroscopic Instrument, AJ 166 (2023) 259 [2306 .06309]. 


J. Guy, S. Bailey, A. Kremin, S. Alam, D.M. Alexander, C. Allende Prieto et al., The 
Spectroscopic Data Processing Pipeline for the Dark Energy Spectroscopic Instrument, AJ 165 
(2023) 144 [2209 . 14482]. 


— 48 — 


105 


106 


107 
108 


109 


110 


111 


112 


113 
114 


[115] 


116 


117 
118 


119 
120 


[121 


[122 


[123 


[124 


A.J. Ross et al., Construction of Large-scale Structure Catalogs for Data from the Dark 
Energy Spectroscopic Instrument, in prep. (2024) [xxxx . xxxxx]. 


A. Rosado-Marin et al., Mitigating Imaging Systematics for DESI DR1 Emission Line 
Galaxies and Beyond, in preparation (2024) . 


E. Chaussidon et al., DESI Y1 fNL, in prep. (2024) [xxxx . xxxxx]. 


J. Yu et al., ELG Spectroscopic Systematics Analysis of the DESI Data Release 1, in 
preparation (2024) . 


A. Krolewski et al., /mpact and mitigation of spectroscopic systematics on DESI Y1 
clustering, in preparation (2024) . 


J. Lasker, A.C. Rosell, A.D. Myers, A.J. Ross, D. Bianchi, M.M.S. Hanif et al., Production of 
Alternate Realizations of DESI Fiber Assignment for Unbiased Clustering Measurement in 
Data and Simulations, arXiv e-prints (2024) arXiv:2404.03006 [2404 . 03006]. 


D. Bianchi et al., Characterization of DESI fiber assigmment incompleteness effect on 2-point 
clustering and mitigation methods for 2024 analysis, in preparation (2024) . 


M. Pinon et al., Mitigation of DESI fiber assigmment incompleteness effect on two-point 
clustering with small angular scale truncated estimators, in preparation (2024) . 


C. Zhao et al., Simulations of DESI Y1 , in prep. (2024) [xxxx.xxxxx]. 


A. Bault, D. Kirkby, J. Guy, A. Brodzeller, J. Aguilar, S. Ahlen et al., Impact of Systematic 
Redshift Errors on the Cross-correlation of the Lyman-a Forest with Quasars at Small Scales 
Using DESI Early Data, arXiv e-prints (2024) arXiv:2402.18009 [2402. 18009]. 


C. Ramirez-Pérez, I. Pérez-Rafols, A. Font-Ribera, M.A. Karim, E. Armengaud, J. Bautista 
et al., The Lyman-a forest catalogue from the Dark Energy Spectroscopic Instrument Early 
Data Release, MNRAS 528 (2024) 6666 [2306 . 06312]. 


S. Filbert, P. Martini, K. Seebaluck, L. Ennesser, D.M. Alexander, A. Bault et al., Broad 
Absorption Line Quasars in the Dark Energy Spectroscopic Instrument Early Data Release, 
arXiv e-prints (2023) arXiv:2309.03434 [2309 . 03434]. 


J. Zou et al., The DESI Damped Lya System Survey: Data Release 1, in preparation (2024) . 


Gaia Collaboration, T. Prusti, J.H.J. de Bruijne, A.G.A. Brown, A. Vallenari, C. Babusiaux 
et al., The Gaia mission, A&A 595 (2016) Al [1609.04153]. 


J. Moustakas et al., fastspecfit , in prep. (2024) [xxxx.xxxxx]. 


J. Moustakas, D. Scholte, B. Dey and A. Khederlarian, “FastSpecFit: Fast spectral synthesis 
and emission-line fitting of DESI spectra.” Astrophysics Source Code Library, record 
ascl:2308.005, Aug., 2023. 


E.L. Wright, P.R.M. Eisenhardt, A.K. Mainzer, M.E. Ressler, R.M. Cutri, T. Jarrett et al., 
The Wide-field Infrared Survey Explorer (WISE): Mission Description and Initial On-orbit 
Performance, AJ 140 (2010) 1868 [1008.0031]. 


D. Valcin et al., Combined tracer analysis for DESI 2024 BAO analysis, in preparation (2024) 


A. Heinesen, C. Blake and D.L. Wiltshire, Quantifying the accuracy of the Alcock-Paczyríski 
scaling of baryon acoustic oscillation measurements, J. Cosmology Astropart. Phys. 2020 
(2020) 038 [1908.11508]. 


J.L. Bernal, T.L. Smith, K.K. Boddy and M. Kamionkowski, Robustness of baryon acoustic 
oscillation constraints for early- Universe modifications of ACDM cosmology, Phys. Rev. D 
102 (2020) 123515 [2004 . 07263]. 


— 49 — 


125 


126 


127 


128 


[129 


130 
131 


132 


133 


[134 


[135 


[136 


[137 


[138 


[139 


[140 


[141 


] 


] 


J. Pan, D. Huterer, F. Andrade-Oliveira and C. Avestruz, Compressed baryon acoustic 
oscillation analysis is robust to modified-gravity models, 2312.05177. 


C. Alcock and B. Paczynski, An evolution free test for non-zero cosmological constant, Nature 
281 (1979) 358. 


N. Padmanabhan and M. White, Constraining anisotropic baryon oscillations, Phys. Rev. D 
77 (2008) 123540 [0804 . 0799]. 


P. Carter, F. Beutler, W.J. Percival, J. DeRose, R.H. Wechsler and C. Zhao, The impact of 
the fiducial cosmology assumption on BAO distance scale measurements, MNRAS 494 (2020) 
2076 [1906 . 03035]. 


A. Lewis, A. Challinor and A. Lasenby, Efficient Computation of Cosmic Microwave 
Background Anisotropies in Closed Friedmann-Robertson- Walker Models, ApJ 538 (2000) 473 
[astro-ph/9911177]. 


D. Eisenstein and M. White, Theoretical uncertainty in baryon oscillations, Phys. Rev. D 70 
(2004) 103523 [astro-ph/0407539]. 


D.J. Eisenstein, H.-J. Seo and M. White, On the Robustness of the Acoustic Scale in the 
Low-Redshift Clustering of Matter, ApJ 664 (2007) 660 [astro-ph/0604361]. 


M. Crocce and R. Scoccimarro, Nonlinear evolution of baryon acoustic oscillations, 
Phys. Rev. D 77 (2008) 023533 [0704.2783]. 


K.T. Mehta, H.-J. Seo, J. Eckel, D.J. Eisenstein, M. Metchnik, P. Pinto et al., Galaxy Bias 
and Its Effects on the Baryon Acoustic Oscillation Measurements, ApJ 734 (2011) 94 
[1104.1178]. 


S.-F. Chen, C. Howlett, M. White, P. McDonald, A.J. Ross, H.-J. Seo et al., Baryon Acoustic 
Oscillation Theory and Modelling Systematics for the DESI 2024 results, arXiv e-prints 
(2024) arXiv:2402.14070 [2402 . 14070]. 


D.J. Eisenstein, H.-J. Seo, E. Sirko and D.N. Spergel, Improving Cosmological Distance 
Measurements by Reconstruction of the Baryon Acoustic Peak, ApJ 664 (2007) 675 
[astro-ph/0604362). 


N. Padmanabhan, X. Xu, D.J. Eisenstein, R. Scalzo, A.J. Cuesta, K.T. Mehta et al., A 2 per 
cent distance to z — 0.35 by reconstructing baryon acoustic oscillations - I. Methods and 
application to the Sloan Digital Sky Survey, MNRAS 427 (2012) 2132 [1202.0090]. 


D.J. Eisenstein, H.-J. Seo, E. Sirko and D.N. Spergel, Improving Cosmological Distance 
Measurements by Reconstruction of the Baryon Acoustic Peak, ApJ 664 (2007) 675 
[astro-ph/0604362). 


N. Padmanabhan, X. Xu, D.J. Eisenstein, R. Scalzo, A.J. Cuesta, K.T. Mehta et al., A 2 per 
cent distance to z — 0.35 by reconstructing baryon acoustic oscillations - I. Methods and 
application to the Sloan Digital Sky Survey, MNRAS 427 (2012) 2132 [1202.0090]. 


S. Alam, M. Aubert, S. Avila, C. Balland, J.E. Bautista, M.A. Bershady et al., Completed 
SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological implications from 
two decades of spectroscopic surveys at the Apache Point Observatory, Phys. Rev. D 103 
(2021) 083533 [2007 .08991]. 


E. Paillas, Z. Ding, X. Chen, H. Seo, N. Padmanabhan, A. de Mattia et al., Optimal 
Reconstruction of Baryon Acoustic Oscillations for DESI 2024, arXiv e-prints (2024) 
arXiv:2404.03005 [2404 . 03005]. 


M. Rashkovetskyi, D. Forero-Sánchez, A. de Mattia, D.J. Eisenstein, N. Padmanabhan, 
H. Seo et al., Semi-analytical covariance matrices for two-point correlation function for DESI 
2024 data, arXiv e-prints (2024) arXiv:2404.03007 [2404 . 03007]. 


— 50 — 


[142] 


[143 


[144 


145 


146 


147 


[148] 


149 


150 


151 


152 


153 


154 


155 


156 


157 


158 


159 


O. Alves et al., Analytical covariance matrices of DESI galaxy power spectra, in preparation 
(2024) . 


D. Forero-Sanchez et al., Analytical and EZmock covariance validation for the DESI 2024 
results, in preparation (2024) . 


J. Mena-Fernández, C. Garcia-Quintero, S. Yuan, B. Hadzhiyska, O. Alves, M. Rashkovetskyi 
et al., HOD-Dependent Systematics for Luminous Red Galaxies in the DESI 2024 BAO 
Analysis, arXiv e-prints (2024) arXiv:2404.03008 [2404.03008]. 


C. Garcia-Quintero, J. Mena-Fernández, A. Rocher, S. Yuan, B. Hadzhiyska, O. Alves et al., 
HOD-Dependent Systematics in Emission Line Galaxies for the DESI 2024 BAO analysis, 
arXiv e-prints (2024) arXiv:2404.03009 [2404 . 03009]. 


C. Gordon, A. Cuceu, J. Chaves-Montero, A. Font-Ribera, A.X. González-Morales, J. Aguilar 
et al., 3D correlations in the Lyman-a forest from early DESI data, J. Cosmology Astropart. 
Phys. 2023 (2023) 045 [2308. 10950]. 


A. Cuceu, H.K. Herrera-Alcantar, C. Gordon, P. Martini, J. Guy, A. Font-Ribera et al., 
Validation of the DESI 2024 Lya forest BAO analysis using synthetic datasets, arXiv e-prints 
(2024) arXiv:2404.03004 [2404 . 03004]. 


H.K. Herrera-Alcantar, A. Muñoz-Gutiérrez, T. Tan, A.X. González-Morales, A. Font-Ribera, 
J. Guy et al., Synthetic spectra for Lyman-a forest analysis in the Dark Energy Spectroscopic 
Instrument, arXiv e-prints (2023) arXiv:2401.00303 [2401 . 00303]. 


J. Guy, S.G.A. Gontcho, E. Armengaud, A. Brodzeller, A. Cuceu, A. Font-Ribera et al., 
Characterization of contaminants in the Lyman-alpha forest auto-correlation with DESI, 
arXiv e-prints (2024) arXiv:2404.03003 [2404 . 03003]. 


S. Brieden, H. Gil-Marín, L. Verde and J.L. Bernal, Blind Observers of the Sky, J. Cosmology 
Astropart. Phys. 2020 (2020) 052 [2006 . 10857]. 


U. Andrade, J. Mena-Fernández, H. Awan, A.J. Ross, S. Brieden, J. Pan et al., Validating the 
Galaxy and Quasar Catalog-Level Blinding Scheme for the DESI 2024 analysis, arXiv e-prints 
(2024) arXiv:2404.07282 [2404 . 07282]. 


R.J. Cooke, M. Pettini and C.C. Steidel, One Percent Determination of the Primordial 
Deuterium Abundance, ApJ 855 (2018) 102 [1710.11129]. 


E. Aver, K.A. Olive and E.D. Skillman, The effects of He I 410830 on helium abundance 
determinations, J. Cosmology Astropart. Phys. 2015 (2015) 011 [1503.08146]. 


E. Aver, D.A. Berg, A.S. Hirschauer, K.A. Olive, R.W. Pogge, N.S.J. Rogers et al., A 
comprehensive chemical abundance analysis of the extremely metal poor Leoncino Dwarf 


galaxy (AGC 198691), MNRAS 510 (2022) 373 [2109 00178]. 


R.L. Workman, V.D. Burkert, V. Crede, E. Klempt, U. Thoma, L. Tiator et al., Review of 
Particle Physics, Progress of Theoretical and Experimental Physics 2022 (2022) 083C01. 


N. Schóneberg, The 2024 BBN baryon abundance update, arXiv e-prints (2024) 
arXiv:2401.15054 [2401 . 15054]. 


A.-K. Burns, T.M.P. Tait and M. Valli, PRyMordial: the first three minutes, within and 
beyond the standard model, European Physical Journal C 84 (2024) 86 [2307 . 07061]. 


V. Mossa, K. Stóckel, F. Cavanna, F. Ferraro, M. Aliotta, F. Barile et al., The baryon density 
of the Universe from an improved rate of deuterium burning, Nature 587 (2020) 210. 


Planck Collaboration, N. Aghanim, Y. Akrami, F. Arroja, M. Ashdown, J. Aumont et al., 
Planck 2018 results. I. Overview and the cosmological legacy of Planck, A&A 641 (2020) Al 
[1807 .06205]. 


—51- 


[160] Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi 
et al., Planck 2018 results. V. CMB power spectra and likelihoods, A&A 641 (2020) A5 
[1907 . 12875]. 


Planck Collaboration, Y. Akrami, K.J. Andersen, M. Ashdown, C. Baccigalupi, M. Ballardini 
et al., Planck intermediate results. LVII. Joint Planck LFI and HFI data processing, A&A 
643 (2020) A42 [2007 .04997]. 


G. Efstathiou and S. Gratton, A Detailed Description of the CAMSPEC Likelihood Pipeline 
and a Reanalysis of the Planck High Frequency Maps, The Open Journal of Astrophysics 4 
(2021) 8. 


[161 


[162 


[163 


E. Rosenberg, S. Gratton and G. Efstathiou, CMB power spectra and cosmological parameters 
from Planck PR4 with CamSpec, MNRAS 517 (2022) 4620 [220510869]. 


M. Tristram, A.J. Banday, M. Douspis, X. Garrido, K.M. Górski, S. Henrot-Versillé et al., 
Cosmological parameters derived from the final (PR4) Planck data release, arXiv e-prints 
(2023) arXiv:2309.10034 [2309 . 10034]. 


[164 


[165 


J. Torrado and A. Lewis, Cobaya: code for Bayesian analysis of hierarchical physical models, 
J. Cosmology Astropart. Phys. 2021 (2021) 057 [2005.05290]. 


[166 


J. Torrado and A. Lewis, ^Cobaya: Bayesian analysis in cosmology." Astrophysics Source 
Code Library, record ascl:1910.019, Oct., 2019. 


Planck Collaboration, P.A.R. Ade, N. Aghanim, C. Armitage-Caplan, M. Arnaud, 
M. Ashdown et al., Planck 2013 results. XVII. Gravitational lensing by large-scale structure, 
A&A 571 (2014) A17 [1303.5077]. 


Planck Collaboration, P.A.R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont et al., 
Planck 2015 results. XV. Gravitational lensing, A&A 594 (2016) A15 [1502.01591]. 


[167 


[168 


[169 


J. Carron, M. Mirmelstein and A. Lewis, CMB lensing from Planck PR4 maps, J. Cosmology 
Astropart. Phys. 2022 (2022) 039 [2206 .07773]. 


F.J. Qu, B.D. Sherwin, M.S. Madhavacheril, D. Han, K.T. Crowley, I. Abril-Cabezas et al., 
The Atacama Cosmology Telescope: A Measurement of the DR6 CMB Lensing Power 
Spectrum and its Implications for Structure Growth, arXiv e-prints (2023) arXiv:2304.05202 
[2304 . 05202]. 


171] N. MacCrann, B.D. Sherwin, F.J. Qu, T. Namikawa, M.S. Madhavacheril, I. Abril-Cabezas 
et al., The Atacama Cosmology Telescope: Mitigating the impact of extragalactic foregrounds 
for the DR6 CMB lensing analysis, arXiv e-prints (2023) arXiv:2304.05196 [2304 . 05196]. 


172] G. Efstathiou, W.J. Sutherland and S.J. Maddox, The cosmological constant and cold dark 
matter, Nature 348 (1990) 705. 


173] S.D.M. White, J.F. Navarro, A.E. Evrard and C.S. Frenk, The Baryon content of galaxy 
clusters: A Challenge to cosmological orthodoxy, Nature 366 (1993) 429. 


[170 


174] D. Scolnic et al., The Pantheon+ Analysis: The Full Data Set and Light-curve Release, 
Astrophys. J. 938 (2022) 113 [2112.03863]. 


175] E.R. Peterson et al., The Pantheon+ Analysis: Evaluating Peculiar Velocity Corrections in 
Cosmological Analyses with Nearby Type la Supernovae, Astrophys. J. 938 (2022) 112 
[2110.03487]. 


[176 


A. Lewis, A. Challinor and A. Lasenby, Efficient Computation of Cosmic Microwave 
Background Anisotropies in Closed Friedmann-Robertson- Walker Models, ApJ 538 (2000) 473 
[astro-ph/9911177]. 


[177] C. Howlett, A. Lewis, A. Hall and A. Challinor, CMB power spectrum parameter degeneracies 


— 52 — 


178 


179 


180 


181 


182 


183 


184 
185 
186 


187 


[188 


[189 


[190 


[191 


[192 


[193 


[194 


[195 


[196 


] 


] 


in the era of precision cosmology, J. Cosmology Astropart. Phys. 2012 (2012) 027 
[1201 . 3654]. 


A. Lewis and S. Bridle, Cosmological parameters from CMB and other data: A Monte Carlo 
approach, Phys. Rev. D 66 (2002) 103511 [astro-ph/0205436]. 


A. Lewis, Efficient sampling of fast and slow cosmological parameters, Phys. Rev. D 87 (2013) 
103529 [1304. 4473]. 


M. Chevallier and D. Polarski, Accelerating Universes with Scaling Dark Matter, International 
Journal of Modern Physics D 10 (2001) 213 [gr-qc/0009008]. 


E.V. Linder, Exploring the Expansion History of the Universe, Phys. Rev. Lett. 90 (2003) 
091301 [astro-ph/0208512]. 


R.M. Neal, Taking Bigger Metropolis Steps by Dragging Fast Variables, arXiv Mathematics 
e-prints (2005) math/0502099 [math/0502099]. 


A. Gelman and D.B. Rubin, Inference from Iterative Simulation Using Multiple Sequences, 
Statistical Science 7 (1992) 457 . 


A. Lewis, GetDist: a Python package for analysing Monte Carlo samples, 1910.13970. 
H. Dembinski and P.O. et al., scikit-hep/iminuit, . 


F. James and M. Roos, Minuit: A System for Function Minimization and Analysis of the 
Parameter Errors and Correlations, Comput. Phys. Commun. 10 (1975) 343. 


A. Cuceu, J. Farr, P. Lemos and A. Font-Ribera, Baryon Acoustic Oscillations and the Hubble 
constant: past, present and future, J. Cosmology Astropart. Phys. 2019 (2019) 044 
[1906. 11628]. 


K.S. Dawson, D.J. Schlegel, C.P. Ahn, S.F. Anderson, É. Aubourg, S. Bailey et al., The 
Baryon Oscillation Spectroscopic Survey of SDSS-III, AJ 145 (2013) 10 [1208.0022]. 


K.S. Dawson, J.-P. Kneib, W.J. Percival, S. Alam, F.D. Albareti, S.F. Anderson et al., The 
SDSS-IV Extended Baryon Oscillation Spectroscopic Survey: Overview and Early Data, AJ 
151 (2016) 44 [1508 . 04473]. 


D.G. York, J. Adelman, J. Anderson, John E., S.F. Anderson, J. Annis, N.A. Bahcall et al., 
The Sloan Digital Sky Survey: Technical Summary, AJ 120 (2000) 1579 [astro-ph/0006396]. 


F. Beutler, C. Blake, M. Colless, D.H. Jones, L. Staveley-Smith, L. Campbell et al., The 6dF 
Galaxy Survey: baryon acoustic oscillations and the local Hubble constant, MNRAS 416 
(2011) 3017 [1106.3366]. 


E.A. Kazin, J. Koda, C. Blake, N. Padmanabhan, S. Brough, M. Colless et al., The WiggleZ 
Dark Energy Survey: improved distance measurements to z — 1 with reconstruction of the 
baryonic acoustic feature, MNRAS 441 (2014) 3524 [1401 .0358]. 


DES COLLABORATION collaboration, Dark Energy Survey Year 1 Results: A Precise H0 
Estimate from DES Y1, BAO, and D/H Data, Mon. Not. Roy. Astron. Soc. 480 (2018) 3879 
[1711.00403]. 


DES collaboration, Assessing tension metrics with dark energy survey and Planck data, Mon. 
Not. Roy. Astron. Soc. 505 (2021) 6179 [2012.09554]. 


W. Sutherland, On measuring the absolute scale of baryon acoustic oscillations, Mon. Not. 
Roy. Astron. Soc. 426 (2012) 1280 [1205.0715]. 


BOSS COLLABORATION collaboration, Cosmological implications of baryon acoustic 
oscillation measurements, Phys. Rev. D 92 (2015) 123516. 


— 53 - 


197 


198 


199 


200 


201 


202 


203 


204 


205 


206 


207 


208 


209 
210 


211 


212 


213 


214 


215 


216 
217 


218 


A.J. Cuesta, L. Verde, A. Riess and R. Jimenez, Calibrating the cosmic distance scale ladder: 
the role of the sound-horizon scale and the local expansion rate as distance anchors, MNRAS 
448 (2015) 3463 [1411.1094]. 


G.E. Addison, G. Hinshaw and M. Halpern, Cosmological constraints from baryon acoustic 
oscillations and clustering of large-scale structure, MNRAS 436 (2013) 1674 [1304.6984]. 


S. Nadathur, W.J. Percival, F. Beutler and H.A. Winther, Testing Low-Redshift Cosmic 
Acceleration with Large-Scale Structure, Phys. Rev. Lett. 124 (2020) 221301 [2001 . 11044]. 


E. Di Valentino, A. Melchiorri and J. Silk, Planck evidence for a closed Universe and a 
possible crisis for cosmology, Nature Astronomy 4 (2020) 196 [1911 .02087]. 


W. Handley, Curvature tension: Evidence for a closed universe, Phys. Rev. D 103 (2021) 
L041301 [1908.09139]. 


J. Martin, Everything you always wanted to know about the cosmological constant problem 
(but were afraid to ask), Comptes Rendus Physique 13 (2012) 566 [1205.3365]. 


C.P. Burgess, The Cosmological Constant Problem: Why it's hard to get Dark Energy from 
Micro-physics, arXiv e-prints (2013) arXiv:1309.4133 [1309.4133]. 


A. Padilla, Lectures on the Cosmological Constant Problem, arXiv e-prints (2015) 
arXiv:1502.05296 [1502.05296]. 


E.V. Linder, The Mirage of w=-1, arXiv e-prints (2007) arXiv:0708.0024 [0708 . 0024]. 


J.A. Frieman, D. Huterer, E.V. Linder and M.S. Turner, Probing dark energy with 
supernovae: Exploiting complementarity with the cosmic microwave background, Phys. Rev. D 
67 (2003) 083505 [astro-ph/0208100]. 


R. de Putter and E.V. Linder, Calibrating dark energy, J. Cosmology Astropart. Phys. 2008 
(2008) 042 [0808.0189]. 


W. Fang, W. Hu and A. Lewis, Crossing the phantom divide with parametrized 
post-Friedmann dark energy, Phys. Rev. D 78 (2008) 087303 [0808.3125]. 


H. Jeffreys, The Theory of Probability, Oxford Classic Texts in the Physical Sciences (1939). 


R. Trotta, Applications of Bayesian model selection to cosmological parameters, Mon. Not. 
Roy. Astron. Soc. 378 (2007) 72 [astro-ph/0504022]. 


W.J. Handley, M.P. Hobson and A.N. Lasenby, PolyChord: nested sampling for cosmology, 
Mon. Not. Roy. Astron. Soc. 450 (2015) L61 [1502.01856]. 


W. Handley, anesthetic: nested sampling visualisation, J. Open Source Softw. 4 (2019) 1414 
[1905 . 04768]. 


A.R. Liddle, Information criteria for astrophysical model selection, Mon. Not. Roy. Astron. 
Soc. 377 (2007) L74 [astro-ph/0701113]. 


S. Grandis, D. Rapetti, A. Saro, J.J. Mohr and J.P. Dietrich, Quantifying tensions between 
CMB and distance data sets in models with free curvature or lensing amplitude, Mon. Not. 
Roy. Astron. Soc. 463 (2016) 1416 [1604.06463]. 


D. Huterer and M.S. Turner, Probing the dark energy: Methods and strategies, Phys. Rev. D 
64 (2001) 123527 [astro-ph/0012510]. 


A. Albrecht et al., Report of the Dark Energy Task Force, astro-ph/0609591. 


E. Hubble, A Relation between Distance and Radial Velocity among Extra- Galactic Nebulae, 
Proceedings of the National Academy of Science 15 (1929) 168. 


G. de Vaucouleurs and G. Bollinger, The extragalactic distance scale. VII - The 
velocity-distance relations in different directions and the Hubble ratio within and without the 
local supercluster, ApJ 233 (1979) 433. 


— 54 — 


219 


220 


221 


222 


[223 


[224] 


[225] 


226 
227 


228 


229 
230 


231 


232 


233 


234 


235 


[236 


A. Sandage and G.A. Tammann, Steps toward the Hubble constant. VIII - The global value, 
ApJ 256 (1982) 339. 


W.L. Freedman, B.F. Madore, B.K. Gibson, L. Ferrarese, D.D. Kelson, S. Sakai et al., Final 
Results from the Hubble Space Telescope Key Project to Measure the Hubble Constant, ApJ 
553 (2001) 47 [astro-ph/0012376]. 


Planck Collaboration, P.A.R. Ade, N. Aghanim, C. Armitage-Caplan, M. Arnaud, 
M. Ashdown et al., Planck 2013 results. XVI. Cosmological parameters, A&A 571 (2014) A16 
[1303 . 5076]. 


A.G. Riess, L.M. Macri, S.L. Hoffmann, D. Scolnic, S. Casertano, A.V. Filippenko et al., A 
2.4% Determination of the Local Value of the Hubble Constant, ApJ 826 (2016) 56 
[1604 .01424]. 


A.G. Riess, S. Casertano, W. Yuan, L. Macri, J. Anderson, J.W. MacKenty et al., New 
Parallaxes of Galactic Cepheids from Spatially Scanning the Hubble Space Telescope: 
Implications for the Hubble Constant, ApJ 855 (2018) 136 [1801.01120]. 


D.W. Pesce, J.A. Braatz, M.J. Reid, A.G. Riess, D. Scolnic, J.J. Condon et al., The 
Megamaser Cosmology Project. XIII. Combined Hubble Constant Constraints, ApJ 891 (2020) 
L1 [2001 .09213]. 


A.G. Riess et al., A Comprehensive Measurement of the Local Value of the Hubble Constant 
with 1 km s~! Mpc-! Uncertainty from the Hubble Space Telescope and the SHOES Team, 
Astrophys. J. Lett. 934 (2022) L7 [2112.04510]. 


W.L. Freedman, B.F. Madore, T. Hoyt, LS. Jang, R. Beaton, M.G. Lee et al., Calibration of 
the Tip of the Red Giant Branch, ApJ 891 (2020) 57 [2002.01550]. 


L. Verde, T. Treu and A.G. Riess, Tensions between the early and late Universe, Nature 
Astronomy 3 (2019) 891 [1907 . 10625]. 


E. Di Valentino, O. Mena, S. Pan, L. Visinelli, W. Yang, A. Melchiorri et al., /n the realm of 
the Hubble tension-a review of solutions, Classical and Quantum Gravity 38 (2021) 153001 
[2103.01183]. 


M. Kamionkowski and A.G. Riess, The Hubble Tension and Early Dark Energy, Annual 
Review of Nuclear and Particle Science 73 (2023) 153 [2211.04492]. 


D.J. Fixsen, The Temperature of the Cosmic Microwave Background, ApJ 707 (2009) 916 
[0911.1955]. 


KATRIN Collaboration, The KATRIN experiment - a direct measurement of the electron 
antineutrino mass in the sub-eV region, Nucl. Phys. A 752 (2005) 14. 


M. Aker, A. Beglarian, J. Behrens, A. Berlev, U. Besserer, B. Bieringer et al., First direct 
neutrino-mass measurement with sub-eV sensitivity, arXiv e-prints (2021) arXiv:2105.08533 
[2105 . 08533]. 


M.C. Gonzalez-Garcia, M. Maltoni and T. Schwetz, NuFIT: Three-Flavour Global Analyses of 
Neutrino Oscillation Experiments, Universe 7 (2021) 459 [2111 .03086]. 


J. Lesgourgues and S. Pastor, Massive neutrinos and cosmology, Phys. Rep. 429 (2006) 307 
[astro-ph/0603494]. 


E. Di Valentino, T. Brinckmann, M. Gerbino, V. Poulin, F.R. Bouchet, J. Lesgourgues et al., 
Exploring cosmic origins with CORE: Cosmological parameters, J. Cosmology Astropart. 
Phys. 2018 (2018) 017 [1612.00021]. 


S.R. Choudhury and $. Hannestad, Updated results om neutrino mass and mass hierarchy 
from cosmology with Planck 2018 likelihoods, J. Cosmology Astropart. Phys. 2020 (2020) 037 
[1907 . 12598]. 


—55-— 


237| DESI Collaboration, DESI 2024 VII: Cosmological constraints from full-shape analyses of the 
two-point clustering statistics measurements, in preparation (2024) . 


238] N. Palanque-Delabrouille, C. Yeche, N. Schóneberg, J. Lesgourgues, M. Walther, S. Chabanier 
et al., Hints, neutrino bounds, and WDM constraints from SDSS DR14 Lyman-a and Planck 
full-survey data, J. Cosmology Astropart. Phys. 2020 (2020) 038 [1911.09073]. 


239] E. Di Valentino, S. Gariazzo and O. Mena, Most constraining cosmological neutrino mass 
bounds, Phys. Rev. D 104 (2021) 083504 [2106 . 15267]. 


240| S. Brieden, H. Gil-Marín and L. Verde, Model-agnostic interpretation of 10 billion years of 
cosmic evolution traced by BOSS and eBOSS data, J. Cosmology Astropart. Phys. 2022 
(2022) 024 [2204 . 11868]. 


241] R.L. Workman, V.D. Burkert, V. Crede, E. Klempt, U. Thoma, L. Tiator et al., Review of 
Particle Physics, Progress of Theoretical and Experimental Physics 2022 (2022) 083C01. 


242] A.F. Heavens and E. Sellentin, Objective Bayesian analysis of neutrino masses and hierarchy, 
J. Cosmology Astropart. Phys. 2018 (2018) 047 [1802.09450]. 


243| R. Jimenez, C. Pena-Garay, K. Short, F. Simpson and L. Verde, Neutrino masses and mass 
hierarchy: evidence for the normal hierarchy, J. Cosmology Astropart. Phys. 2022 (2022) 006 
[2203 . 14247]. 


244| S. Gariazzo, M. Gerbino, T. Brinckmann, M. Lattanzi, O. Mena, T. Schwetz et al., Neutrino 
mass and mass ordering: no conclusive evidence for normal ordering, J. Cosmology 
Astropart. Phys. 2022 (2022) 010 [2205 . 02195]. 


245] D. Baumann, D. Green and M. Zaldarriaga, Phases of New Physics in the BAO Spectrum, J. 
Cosmology Astropart. Phys. 2017 (2017) 007 [1703.00894]. 


246] D. Baumann, F. Beutler, R. Flauger, D. Green, A. Slosar, M. Vargas-Magaña et al., First 
constraint on the neutrino-induced phase shift in the spectrum of baryon acoustic oscillations, 
Nature Physics 15 (2019) 465 [1803.10741]. 


247] A. Perez-Fernandez, R. Ruggeri et al., Fiducial Cosmology systematics for DESI 2024 BAO 
Analysis, in preparation (2024) . 


248] DESI Collaboration, DESI 2024 VIII: Constraints on Primordial Non-Gaussianities, in 
preparation (2024) . 


A Comparison of results from DESI and combined DESI--SDSS BAO data 


As discussed in Section 3.3, it is possible to define a combined “DESI+SDSS” BAO dataset 
across the redshift range covered by both surveys which uses: 1) SDSS results at zef = 
0.15,0.38 and 0.51 in place of the DESI BGS and lowest-redshift LRG points; 2) DESI 
results from LRGs in 0.6 < z < 0.8, the combination of LRGs and ELGs in 0.8 < z < 1.1, 
and ELGs and QSOs at higher redshifts; and 3) the combined DESI+SDSS result from 
Eqs. (3.3) and (3.4) for the Lya BAO. The combined DESI+SDSS dataset therefore leverages 
observations from SDSS and DESI in the redshift regimes where each survey respectively 
covers a larger effective volume, and combines the Lya BAO measurements from both. We 
emphasise again that unlike the DESI DR1 BAO, this combination of measurements does 
not come from uniform analysis methods or data processing pipelines; nevertheless, they are 
of interest as, by maximising the effective volume at each redshift, they provide the most 
precise combination of BAO data currently available. In this section we compare the results 
using this combination to the DESI ones in the main body of the paper, in the cases of most 
interest. 
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We start with the flat ACDM model. DESI+SDSS constraints alone give 


Qm = 0.297 + 0.012, 
rah = (101.3 + 1.1) Mpc, 


| (DESI+SDSS). (A.1) 


This is in excellent agreement with the DESI-only results in Eq. (4.1), with about 20% smaller 
errors. Combining BAO information with BBN, we obtain 


Ho = (67.98 + 0.75)kms~'Mpe~! — ((DESL-SDSS) BAO+BBN), (A.2) 


and further adding the acoustic angular scale constraint, 


Ho —(68.132:0.59) kms ! Mpe !  ((DESI+SDSS)+BBN+9,), (A.3) 


both of which are in excellent agreement with the corresponding DESI results (Eqs. (4.4) 
and (4.5)), and with those from the CMB. These results are also shown in Figure 9. The 
combination with the full CMB data gives 


Qm = 0.309 + 0.0048, 


Ho = (67.80 + 0.37) kms"! Mpc”! | (DESI+SDSS) BAO+CMB. (A-4) 


Finally, the 95% upper limit on the neutrino mass in flat ACDM for the combination of this 
BAO data with the CMB, and with the >) m, > 0 eV prior is 


Y m, «0.082eV (95%, (DESI+SDSS) BAO+CMB). (A.5) 


This is slightly weaker than the upper bound in Eq. (7.2): as discussed in Section 7.1, it is 
because the BAO precision in both cases is close enough that the primary determinant of 
the limit obtained on * 5 m, is actually the central value of Ho, a lower Ho allowing a larger 
>) mu. As can be seen from the results quoted above, the DESI+SDSS combination favours 
very slightly lower values of Ho. 

In terms of extensions to the base ACDM model, the case of most interest is naturally 
the (wo, Wa) parametrization of the dark energy equation of state parameter. Figure 12 shows 
the marginalized posterior constraints in wo and Wa for the same combinations with CMB 
and SN Ia datasets as in the right panel of Figure 6, but replacing the DESI DR1 BAO data 
with the DESI+SDSS equivalents. We find 


wo = 0.855 + 0.000, | (DESI+SDSS)+CMB (A.6) 
Wa = —0.607) 38, +PantheonPlus, 

from combination with PantheonPlus, 
wo = —0.692 + 0.095, | (DESI+SDSS)+CMB (A.7) 
Wa = —1.06*038, +Union3, 


= 5? = 


(DESI+SDSS) BAO + CMB + PantheonPlus 


= 
E (DESI+SDSS) BAO + CMB + Union3 
EE (DESI+SDSS) BAO + CMB + DESY5 


—1.0 —0.8 —0.6 —0.4 

wo 
Figure 12. The 68% and 95% marginalized posterior constraints from DESI+SDSS BAO combined 
with CMB and each of the PantheonPlus, Union3 and DESY5 SN Ia datasets. Compared to cases 


with DESI DR1 BAO shown in the right panel of Figure 6 all contours are very consistent but the 
significances of the tensions with ACDM marginally decrease. 


with Union3 SN Ia, and 


Wa = —0.88* 022. +DESY5, 


when using the DESY5 SN la data, all in excellent agreement with the DESI results provided 
in the main text. All these results shift marginally closer to the ACDM expectation: the 
AX? jap values between the maximum a posteriori of the wowaCDM model and the maximum 
of the posterior fixing (wo, wa) = (—1,0) are —7.0, —12.6 and —15.5 for the combinations 
with PantheonPlus, Union3 and DESY5 respectively, corresponding to discrepancies at the 
2.20, 3.10 and 3.50 significance levels. These are all comparable to, though slightly lower 
than, the equivalent 2.5c, 3.50 and 3.90 discrepancies reported earlier using DESI DR1 BAO. 

Since replacing all the DESI points at z « 0.6 with SDSS does not change the qualitative 
conclusions of tension with ACDM, this indicates that—despite being a ~ 20 offset from the 
ACDM model—the single DESI BAO measurement of Dy /Dy(z = 0.51) from LRGs is not 
significantly driving the dark energy results in Section 5.2. 

A better understanding of why this is the case can be obtained from Figure 13, which 
shows the DESI BAO measurements of Dy/ra and Dy/Dy relative to the best-fit ACDM 
model, and overlays the model predictions for two representative wow; CDM models. The 
zet = 0.51 DESI BAO measurement is anomalous in the less well-measured quantity Dy / Dg 
(at the S 20 level), although perfectly consistent with ACDM in Dy/rg. The blue dash-dot 
lines show the model which has the minimum Xt Ap for a fit to the DESI BAO alone: this is 
a rather extreme case, with Qm = 0.385, wg = —0.159 and wa = —3.0, the wa value hitting 
our prior Wa € |-3,2] (Table 2). (However, as noted in Section 5.2, although these best-fit 
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—^* wowaCDM fit to DESI alone 
—*: wgw,CDM fit to DESI+CMB+DESY5 


0.5 1.0 1.5 2.0 2.5 0.5 1.0 1.5 2.0 2:5 
Redshift z Redshift z 


Figure 13. Same as the lower panels of Fig. 1, except with model lines for two wow, CDM models 
shown. The blue dot-dashed line shows the best-fit wow; CDM model from DESI data alone, while 
the orange line shows the best fit from DESI+CMB+DESY5. As before, the solid line (at value of 
unity) is the best-fit ACDM model to DESI data alone, while the dashed line and the grey shaded 
region showing the Planck best-fit ACDM model and the corresponding 68% uncertainties. 


parameters are extreme, the level of *tension" with ACDM actually does not even reach the 
20 level, with Ax?;4p = —3.7.) The orange dash-dot lines instead show the model predictions 
for the model with Qm = 0.316, wo = —0.733 and wa = —1.010, which minimises x?;4p for 
the fit to the combination of DESI BAO+CMB+DESY5. This is the data combination which 
also gives the highest significance tension with ACDM, of 3.90. 

What is clear from Figure 13 is that the slightly anomalous DESI measurement of 
Du/Du(z = 0.51) does not strongly affect the wow; CDM fits: even the (wo, Wa) values at 
the edge of our prior range fail to fit this data point well, so it contributes very little to the 
Ax? between ACDM and wow; CDM even for the fit to DESI BAO alone (in which case the 
statistical evidence for variation of the equation of state is very low anyway). In the more 
interesting case of combining BAO, CMB and SN Ia data the difference between the ACDM 
and woWaCDM model predictions in Dy/Dg is even smaller, due to the tight constraints 
from SN Ia at z « 0.5. 

In conclusion, the DESI result for Dy /Dg(z = 0.51) appears to be simply a mild (~ 20) 
statistical fluctuation. It does not appear to have a simple explanation in terms of any model, 
and it consequently also does not strongly drive the results in any model that we studied, 
especially when DESI data are combined with CMB and SN Ia. We have explicitly shown 
that our conclusions about dark energy still hold when all DESI measurements at z « 0.6 
are completely replaced by the corresponding SDSS values. 


B Comparison to Planck PRA likelihoods 


Throughout the main text of this paper, we have used the plik likelihood from the official 
Planck PR3 data release as the default for the high-Z CMB temperature and polarisation 
(TTTEEE) power spectrum likelihood analysis. As mentioned in Section 2.4.2, more recently 
other likelihoods have been released by teams that make use of updated data from the 
Planck PRA release [162-164]. The use of these likelihoods leads to small shifts in parameter 
estimates compared to plik, but in almost all cases of interest in this paper, after combination 
with CMB lensing likelihoods from Planck and ACT [22, 169], these shifts become negligibly 
small, and even more so when our DESI BAO likelihoods are included. Nevertheless, for 
completeness we report here the results obtained from replacing the plik likelihood with the 
CamSpec likelihood from Refs. [162, 163], restricting attention to only those models where the 
parameter shifts are non-negligible. For conciseness, in this appendix we will refer to shifts 
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in the mean value of a parameter that amount to r% of the original statistical uncertainty as 
a shift of “xo”, but this is not a statement about the level of discrepancy between the results 
(if any), and performing such a calculation accounting for the overlap in the data used by 
the different CMB likelihoods is far beyond the scope of this paper. 

In the base ACDM model, we find 


Ho = 67.73 + 0.36 kms! Mpe”! (DESI BAO+Planck|CamSpec] (B.1) 
+CMB lensing,) 


corresponding to a 0.66 shift in the mean, and marginal improvement in the uncertainties 
compared to the result in Table 3 using plik. Although not quantities we have focused on in 
this paper, the baryon density wp and matter power spectrum amplitude og change slightly 
as well, from 


wp = 0.02248 + 0.00013, DESI BAO+ Planck |plik] (B.2) 
gg = 0.8135 + 0.0053, +CMB lensing, 


to 


Wb = 0.02226 + 0.00012, DESI+ Planck [CamSpec] (B.3) 
gg = 0.8103 + 0.0051, +CMB lensing, 


a shift of ~ 1.80 in wy. This causes a smaller shift in the sound horizon value, from rg = 
147.34 + 0.22 to ra = 147.59 + 0.21 when swapping from plik to CamSpec, although with 
negligible change to rgh. Changes in Qm are also negligible (< 0.50). 

Although it has been noted [163] that the CamSpec likelihood somewhat alleviates the ap- 
parent preference in plik for a non-zero spatial curvature, related to the Ar, lensing amplitude 
parameter anomaly, this preference is in any case also removed by combination with CMB 
lensing and BAO, and we find the result in Eq. (4.10) is completely unaffected by the change 
in Planck likelihood. The upper bounds on 5) m, tighten marginally, e.g., >) m, < 0.071 eV 
in flat ACDM and assuming a prior )) m, > 0 eV, compared to Eq. (7.2). There is also a 
small shift in the result for Neg in flat ACDM, 


Neg = 3.20 + 0.19, (DESI BAO+ Planck |CamSpec]4- (B.4) 
CMB lensing), 


a shift of slightly more than 0.507 and a marginal change in the uncertainty compared to 
Eq. (7.8). 

In models with changes to the dark energy equation of state, we find that changing 
between the plik and CamSpec likelihoods has a negligible difference on wo and wa. 
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